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ABSTRACT 

A  normally  conducting  contrarotating  homopolar  motor  has  been  designed  and 
constructed.  The  reaction  torque,  in  the  outer  rotor,  from  the  inner  rotor  is 
utilized  to  produce  true  contrarotation.  The  machine  utilizes  liquid  cooled 
conductors,  high  performance  liquid  metal  current  collectors,  and  ferrous 
conductors  in  the  active  region.  The  basic  machine  output  is  175  hp  at  ±  1,200 
rpm  with  an  input  of  4  volts  and  35,000  amps. 

ADMINISTRATIVE  INFORMATION 
The  motor  design  effort  described  herein  was  initiated  under  the  Independent 
Research  and  Independent  Exploratory  Development  Program  at  this  center  under  Work 
Unit  1-27 1 1-100  and  subsequently  supported  by  the  Office  of  Naval  Research  under 
Program  Element  62121N  Projects  RH21E41  and  RH21E42.  This  project  milestone 
report  is  submitted  in  partial  fultillment  of  Milestone  1  Task  1  of  the  Advanced  Electrical 
Systems  Project  (RH21E41)  of  the  Surface  Ship  Technology  Program 
(SC1A/PH)602121N).  The  work  described  herein  was  sponsored  by  ONR  334  and  was 
I)erformed  by  personnel  of  the  Carderock  Division,  Naval  Surface  Warfare  Center 
(CDNSWC  811),  Annapolis,  Maryland. 
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INTRODUCTION 

As  a  result  of  numerous  ship  impact  studies  relating  technology  to  performance 
and  other  attributes,  contrarotating  propulsion  trains  have  been  identified  as  viable 
candidates  for  future  surface  combatants  In  addition,  the  homopolar  DC.  motor  has 
been  identified  as  an  attractive  long  term  electric  drive  option  for  both  surface  ships^  and 
submarines^.  Since  the  homopolar  motor  concept  is  well  adapted  to  contrarotation  we 
have  an  opportunity  to  provide  the  advantages  of  both  advanced  electric  drive  and 
contrarotation  at  the  same  time.  In  order  to  establish  the  feasibility  of  this  approach,  a 
development  program  to  demonstrate  direct  drive  contrarotation  on  a  small  scale  was 
undertaken. 
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This  report  documents  the  basic  design,  design  decisions,  and  philosophy  that 
have  gone  into  the  development  and  detailed  design  of  this  175  hp  model  contrarotating 
homopolar  motor.  This  motor  is  being  constructed  as  a  demonstration  both  of 
contrarotation  and  of  normally  conducting  homopolar  motor  technology.  In  addition,  a 
variety  of  design  and  analysis  techniques  have  been  developed  which  will  be  applicable 
to  other  machines,  both  single  and  contrarotating,  and  with  normal  or  superconductive 
excitation. 


BASIC  PRINCIPLES 

This  motor  operates  on  the  same  basic  principles  as  a  single-rotating  homopolar 
motor4.  Torque  is  produced  by  the  direct  interaction  of  current  in  a  conductor  with  a 
separately  produced  magnetic  induction.  The  two  basic  types  of  homopolar  motors  are 
the  disc  machine.  Fig.  1,  and  the  drum  machine.  Fig.  2.  In  a  disc  machine,  conducting 
discs,  with  radial  current,  are  placed  in  an  axial  magnetic  induction.  The  drum  machine 
consists  of  conducting  drums,  with  axial  current,  in  a  radial  magnetic  induction.  Both 
machine  types  are  fitted  with  current  collectors  to  provide  electrical  contact  between  the 
rotating  discs  or  drums  and  the  stationary  current  return  path  that  completes  the  electrical 
circuit  and  allows  connection  to  the  power  source.  This  return  path  most  commonly  takes 
the  same  general  form  as  the  rotor  current  path.  A  disc  machine  has  stationary  discs  that 
make  up  the  return  path  while  a  drum  machine  is  fitted  with  stationary  drums  that  make 
up  the  return  path. 
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Fig.  1.  Single  Turn  Disk  Homopolar  Motor 
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Fig.  2.  Single  Turn  Drum  Homopolar  Motor 
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In  both  types  of  machines  the  current  and  magnetic  induction  interact,  by  means 
of  a  vector  cross  product,  to  produce  a  force  that  is  at  right  angles  to  both  the  current  and 
magnetic  induction  vectors.  Due  to  the  symmetry  of  both  the  magnetic  induction  and  the 
current  carrying  elements  of  the  machine,  the  summation  of  these  forces  results  in  a 
torque  on  the  rotor  elements  and  an  equal  but  opposite  reaction  torque  on  the  stator 
elements  that  make  up  the  stationary  current  return  path. 

In  a  single-rotating  motor,  the  reaction  torque  produced  in  the  stationary  elements 
can  produce  no  power.  In  the  case  of  a  contrarotating  motor.  Fig.  3,  the  normally 
stationary  current  return  path,  now  called  the  outer  rotor,  is  allowed  to  rotate  and  thus  can 
also  produce  power.  The  rotor,  now  called  the  inner  rotor,  continues  to  produce  the  same 
torque  and  power  as  before.  Due  to  the  fact  that  outer  rotor  torque  is  the  reaction  torque 
from  the  inner  rotor,  the  torque  on  the  two  rotors  must  be  equal  in  magnitude  but  opposite 
in  direction.  If  the  applied  loads  are  equal,  then  the  rotors  will  rotate  at  the  same  speed, 
but  in  opposite  directions,  and  will  both  produce  the  same  power,  effectively  doubling  the 
output  power  of  the  motor.  Of  course  this  increased  output  power  is  not  available  for 
free.  The  electrical  input  power  is  also  doubled.  Since  the  differential  speed  between  the 
inner  and  outer  rotors  is  now  twice  as  high,  the  motor  will  produce  twice  the  back  EMF, 
effectively  doubling  the  terminal  voltage  of  the  machine.  Two  additional  current 
collectors,  called  terminal  collectors,  are  required  to  make  electrical  contact  with  the 
outer  rotor.  The  outer  rotor  must  additionally  be  fitted  with  bearings,  seals,  rotating 
instrumentation  hookups,  etc.,  to  accommodate  rotation.  A  stationary  housing  surrounds 
the  outer  rotor  providing  a  mounting  point  for  the  bearing  housings  and  also  functions  as 
a  flux  shield. 
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Fig.  3.  Single  Turn  Contrarotating  Drum  Homopolar  Motor 


DESIGN  APPROACH 

This  motor  is  intended  to  be  a  basic  proof  of  principle  and  demonstration 
machine.  It  is  not  intended  to  demonstrate  full  scale  machine  parameters,  such  as  power 
density  or  rotational  speed,  as  this  would  compromise  the  design  of  a  small  motor, 
resulting  in  excessive  risk  and  increased  cost  The  basic  approach  was  to  design  a 
machine  that  was  large  enough  to  allow  meaningful  evaluation  of  motor  design  and 
operation  yet  small  enough  to  be  constructed  in-house.  In  keeping  with  the  basic  intent 
of  the  machine,  we  also  tried  to  keep  the  motor  simple  to  operate  even  if  it  meant  giving 
up  some  power  density.  One  area  that  was  not  compromised,  however,  is  the  thermal 
design.  The  machine  could  have  been  considerably  simplified  if  it  was  not  fully  cooled, 
since  piost  of  the  current  carrying  conductors  are  rotating  and  thus  require  seals  and  a 
complex  manifold  arrangement  for  coolant.  However  an  un-cooled  motor  would  have 
been  difficult  to  evaluate  since  steady  state  operation  would  be  impossible  to  achieve.  An 
un-cooled  design  also  would  have  been  too  large  a  deviation  from  a  full  scale  machine 
design  approach,  reducing  the  usefulness  of  this  model  in  demonstrating  and  evaluating 
the  required  technologies. 
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MAJOR  DESIGN  DECISIONS  AND  RATIONALE 
MACHINE  TYPE 

The  homopolar  concept  was  chosen  for  this  design  because  of  its  many 
recognized  advantages  as  a  direct  drive  propulsion  motor^>3  and  its  adaptability  to  direct 
contrarotation.  The  relatively  straightforward  magnetic  and  armature  circuits  simplify  the 
outer  rotor.  This  allows  fcH*  the  design  of  a  trae  contrarotating  motor  where  the  out«’ 
rotor  is  driven  by  inner  rotor  reaction  torque. 

Disk 

Disk  type  homopolar  machines  have  been  built  but  they  are  more  difficult  to  cool. 
They  also  present  some  unique  problems  for  liquid  metal  current  collectors  which  must 
operate  at  two  different  radii.  As  a  result  of  the  geometry,  liquid  metal  transfer  between 
collectOTS  has  been  a  problem  requiring  elaborate  recirculation  systems  which  have 
always  been  problematic. 

Drum 

The  drum  type  homopolar  has  been  the  most  successful  design  to  date  for  several 
reasons.  The  liquid  metal  current  collectors  are  all  at  one  radius  and  liquid  metal 
migration  can  easily  be  controlled.  A  subset  of  the  drum  type  machine  is  the  bar  type,  hi 
this  arrangement,  the  drum  set  is  replaced  by  bars  evenly  spaced  around  the  rotor.  This 
construction  is  particularly  good  from  a  cooling  standpoint  Each  bar  can  be  provided 
with  one  or  more  cooling  holes  connected  to  coolant  manifolds  at  each  end  of  the  rotor. 
Electrical  interconnection  is  also  facilitated  with  the  bar  arrangement  Since  the  bars  are 
run  the  entire  length  of  the  rotor,  for  ease  of  cooling,  any  collector  disk  can  be  coimected 
to  the  appropriate  bar.  In  a  motor  with  normally  conducting  excitation  the  bar 
construction  has  the  added  advantage  of  lower  reluctance  in  the  active  region.  This  is 
accomplished  either  by  plying  copper  bars  in  slots  cut  in  a  ferrous  rotor  or  by  using 
ferrous  conductors  in  the  active  region.  For  these  reasons  bar  type  construction  was 
chosen  for  the  model  contrarotating  motor. 

MACHINE  SIZE 

Since  this  is  a  demonstration  machine  and  not  a  scale  model,  the  size  and  power 
are  somewhat  arbitrary.  Therefore  the  size  was  based  on  more  practical  considerations 
rather  than  on  sonte  particular  set  of  parameters  derived  from  a  scale  down.  Based  on 
prior  experience  with  the  DTRC  400  hp  motor^  we  felt  that  a  machine  of  approximately 
the  same  overall  physical  dimensions  would  be  desirable.  This  machine  was  large 
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enough  to  allow  meaningful  evaluation  but  not  so  large  as  to  be  difficult  to  construct  in- 
house.  Rough  calculations  indicated  that  a  machine  with  a  magnetic  circuit  that  produced 
between  0.1  and  0.2  Webers  of  flux  in  the  active  region  would  be  in  the  desired  size 
range.  A  more  detailed  magnetic  design  was  done  for  both  the  0.1  and  0.2  Weber  cases. 
The  0. 1  Weber  design  yielded  a  machine  of  the  desired  physical  size  and  with  reasonable 
electrical  parameters.  The  0.1  Weber  design  was  therefore  selected  for  continuation  to 
detailed  design. 

ASSEMBLY  METHOD 

Thoe  are  two  basic  assembly  options  for  single-rotating  homopolar  machines. 

The  first  consists  of  splitting  the  statOT  clam-shell  fashion.  This  allows  the 
complete  rotor  to  be  placed  into  the  lower  stator.  Split  line  seals  are  applied  to  the  lower 
stator  half.  The  upper  statcn-  half  is  then  placed  over  the  lower  stator-rotor  combination 
and  secured  to  the  lower  stator  half.  Ihis  assembly  is  then  placed  in  the  flux  shield  and 
fitted  with  the  bearings  and  seals. 

The  second  method  is  the  axial  stackup.  Widi  this  assembly  method  the  stator 
need  not  be  split  The  machine  is  generally  assembled  with  the  shaft  oriented  votically. 
The  rotOT  and  stator  components  are  alternately  placed  onto  the  assembly,  generally 
around  the  shaft.  This  assembly  can  then  be  placed  in  the  flux  shield,  returned  to  the 
horizontal,  and  fitted  with  bearings  and  seals. 

A  non  split  assembly  method  was  chosen  for  this  contrarotating  motor  for  several 
reasons.  In  a  contrarotating  motor,  what  would  normally  be  the  stator  in  a  single  rotating 
motor  is  now  the  outer  rotor.  Thus,  it  must  handle  rotational  forces  as  well  as  reaction 
torque,  liiere  was  concern  that  this  combination  of  forces  might  result  in  geometric 
distortion  of  the  outer  rotor.  A  full  scale  machine  would  have  much  lower  centrifugal 
loading  since  it  would  opoate  at  relatively  lower  speeds.  In  addition,  for  a  one  turn 
machine,  this  design  offers  construction  simplicity,  although  assembly  difficulty  is 
increased  somewhat.  The  assembly  /  disassembly  procedure  will  be  somewhat  time 
consuming  and  require  care  and  precision.  Vertical  assembly  fixtures  will  be  required  to 
aid  tile  process. 

EXCITATION 

The  choices  for  excitation  of  a  homopolar  motor  are  superconductive  field  coils  or 
conventional  air  or  liquid  cooled  coils.  Since  this  is  a  demonstration  machine  we  felt  that 
the  simplicity  of  operation  provided  by  a  conventional  liquid  cooled  field  was  more 
important  than  the  higher  power  density  a  superconductive  design  could  provide.  A 
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machine  based  on  a  conventional  field  will  still  address  many  technologies  are 
common  to  both  types  of  contrarotating  homopolars,  ie.  current  collectors,  basic  d<»<;ign 
and  analysis,  construction  and  assembly  techniques.  In  addition,  this  design 
the  unique  aspects  of  non-supeiconducdve  designs,  i.e.  ferrous  flux  path,  rotating  field 
coils,  and  the  use  of  ferrous  conductors  in  the  active  regions. 

GENERAL  MACHINE  DESCRIPTION 
OVERALL  MACHINE  LAYOUT 

All  detail  drawings  are  included  in  appendix  E.  An  assembly  section  of  the 
complete  motor  is  shown  on  drawing  A27-19992-4.  Assembly  sections  of  the  inner  rotor, 
outer  rotor,  and  housing  are  shown  on  drawings  A27-19992-7,  A27-19992-8,  and  A27- 
19992-64  respectively.  Section  views  at  various  axial  positions  are  shown  on  A27- 
19992-2,  A27-19992-3,  A27-19992-5,  and  A27-19992-6.  A27-19992-1  is  a  detaUed 
parts  and  drawing  Ust  showing  the  detail  drawing  for  each  piece.  All  piece  numbers  refer 
to  one  or  more  of  these  drawings  and  are  consistent  among  these  drawings.  Drawings  of 
the  tough  forgings  are  not  included. 

This  machine  is  a  single  turn,  normal  conducting,  contrarotating  homopolar 
motor.  The  motor  will  produce  175  hp  at  ±  1200  r/min.  with  a  nominal  input  of  4  volts 
and  35,000  amperes.  Contrarotation  is  accomplished  by  fitting  the  normally  stationary 
current  return  path  (stator)  widi  beatings,  seals,  and  an  additional  pair  of  current 
collectOTS  which  allow  it  to  rotate  opposite  to  die  inner  rotor.  Magnetic  flux  is  generated 
by  a  pair  of  liquid  cooled  coils,  in  a  quadrapole  configuration,  that  rotate  with  the  inner 
rotor.  Die  motor  has  a  complete  ferrous  path  for  the  magnetic  flux.  Included  in  the 
ferrous  path  are  iron  conductors  in  the  active  regions  of  both  inner  and  outer  rotors. 

Current  transfer  between  inner  and  outer  rotors  and  between  the  outer  rotor  and  terminals 
IS  accomplished  with  close  clearance  liquid  metal  current  collectors.  Both  rotors  are 
fitted  with  rolling  element  bearings  and  high  performance  cover  gas  «‘aig  All 
electrically  active  elements  are  directly  cooled  by  forced  flow  liquid  coolant 

MAGNETIC  CIRCUIT 

The  magnetic  circuit  is  made  up  of  the  two  field  coils,  piece  2,  the  iron  bar  sets  of 
the  inner  and  outer  rotoia,  pieces  8  and  16,  and  the  various  other  ferrous  members  that 
shape  and  direct  the  magnetic  flux  through  the  active  regions  of  the  inner  and  outer 
rotors.  Because  of  the  limited  ampere  turns  available  from  a  conventional  cofl,  ferrous 
material  must  be  provided  around  the  complete  flux  path  to  lower  its  reluctance,  thereby 
increasing  the  flux  generated  by  a  given  number  of  ampere  turns  in  the  field  coils.  The 
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iron  bars  in  the  active  regions  of  both  the  inner  and  outer  rotor  are  a  unique  aq)ect  of  this 
design.  A  more  conventional  approach  in  a  non  superconductive  homopolar  machine  is 
to  place  conducting  bars  in  slots  cut  in  a  solid  fenous  rotor.  The  axial  current  is  then 
carried  by  the  conducting  bars  and  the  radial  flux  is  carried  by  the  iron  between  the  slots. 

ELECTRICAL  CIRCUIT 

The  electrical  circuit  includes  all  those  components  required  to  transport  the  axial 
current  through  the  active  regions  of  both  the  inner  and  outer  rotors.  It  is  comprised  of 
the  two  single  rotating  terminal  current  collectors,  pieces  15, 18, 26,  and  27,  two 
contrarotating  current  collectors,  pieces  9, 12, 14,  and  19,  a  set  of  32  inner  rotor  bars, 
piece  8,  a  set  of  36  outer  rotor  bars,  piece  16.  Both  the  iimer  and  outer  rotor  bar  sets  are 
brazed  to  continuous  copper  rings,  pieces  10  and  1 1,  at  both  ends.  This  braze  joint 
provides  the  electrical  connection  to  the  iron  bars  as  well  as  connecting  die  bars 
electrically  in  parallel.  All  of  these  components  are  connected  electrically  in  series.  It  is 
the  interaction  of  the  axial  current  with  the  radial  magnetic  field  in  the  active  region  that 
produces  the  torque  on  the  inner  rotor  and  the  reaction  torque  on  the  outer  rotor.  The 
basic  electrical  circuit  with  the  current  patii  is  shown  in  Fig.  4. 
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Fig.  4.  Electrical  Circuit  Showing  Current  Path 


COOLING  SYSTEM 

TTic  motor  design  features  liquid  cooling  in  both  the  inner  and  outer  rotors, 
including  the  bars  in  the  active  regions,  the  current  collectors,  and  the  field  coils. 
BaseUne  design  requirements  stipulated  an  armature  current  of  35,000  amperes  at  full 
load,  with  an  overload  requirement  of  100,000  amperes.  In  accordance  with  the  decision 
to  design  the  motor  to  the  overload  conditions  wherever  practicable,  the  thermal  analysis 
was  performed  at  die  100,000  ampere  opiating  condition. 

BEARINGS  AND  COVER  GAS  SEALS 

The  inner  rotor  is  supported  by  rolling  element  bearings,  pieces  46  and  50,  which 
are  rmunted  intenud  to  the  outer  rotor  stub  shafts,  pieces  33  and  35.  TTiese  bearings 
operate  in  a  contrarotating  mode  with  the  outer  and  inner  races  rotating  in  opposite 
directions.  The  outCT  rotor  is  in  turn  supported  by  bearings,  pieces  45  and  49,  in 

carriers,  pieces  24  and  40,  which  are  mounted  in  the  stationary  housings,  pieces  39  and 
These  bearings  are  operated  in  a  single  rotating  mode  with  a  stationary  outer  race, 
Mechatucal  face  seals  with  carbon  faces  runrung  on  AI2O3  coated  runners  are 
required  at  both  ends  of  the  machine  to  seal  the  high  purity  cover  gas  in  tire  motor. 
Contrarotating  seals,  pieces  44  and  48,  are  fitted  between  the  inner  rotor  and  the  outer 
rotor  stub  shafts.  Single  rotating  seals,  pieces  43  and  47,  are  also  fitted  between  the  outer 
rotor  stub  shafts  and  the  bearing  housings.  Self  lubricating  lip  seals,  pieces  103, 104,  and 
105,  are  located  between  the  cover  gas  seals  and  the  bearings  to  insure  that  no  bearing 
lubricant  contaminates  the  dry  gas  seals 

SERVICE  MODULE 

To  reduce  the  risk  that  a  failure  in  any  of  the  ancillary  systems  (c.g.  cooling 
system,  instrumentation,  etc.)  could  cause  damage  to  or  require  disassembly  of  the  motor 
itself,  a  removable  service  module  contains  all  required  field  coil  slip  rings, 
instrumentation  cormections,  and  cooling  fluid  seals  and  urtions.  The  module  is  a  totally 
separate  assembly  bolted  to  the  shafts  and  the  beating  housing.  Detailed  drawings  of  the 
service  module  are  included  in  appendix  F. 

SUP  RINGS 

A  slip  ring  assembly,  to  provide  excitation  current  to  the  rotating  field  coils,  is 
located  on  the  service  end  of  the  machine.  It  is  of  conventional  design  using  copper  slip 
rings  and  metal  graphite  brushes. 
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DETAILED  MACHINE  DESIGN 

MAGNETIC  CIRCUIT 

Normally  the  required  electrical  characteristics  of  the  motor,  along  with  the 
desired  number  of  turns,  determines  the  required  flux.  In  this  case,  the  desired  overall 
physical  size  of  the  motor  dictated  the  desired  flux.  Widi  a  requirement  for  0.1  Webo^ 
total  flux  we  could  proceed  with  the  magnetic  design.  Since  this  is  a  single  turn  machine, 
a  quadrapole  configuration  was  chosen,  as  it  provides  a  symmetric  geometry.  Each  half 
of  the  magnetic  circuit  needed  to  produce  0.05  Webers  of  flux.  Using  the  required  flux 
and  a  flux  density  of  1.8  Tesla,  a  magnetic  circuit  Aat  provided  a  constant  cross  sectional 
area  was  designed.  Hie  required  ariqiere  turns  in  the  coil  was  calculated  using  the  total 
path  length  in  ircm  and  in  air.  This  tough  design  was  modeled  in  two  dimensions  with 
finite  element  software.  Based  on  the  finite  element  results,  the  coil  size  and  iron 
geometry  were  adjusted  to  provide  a  relatively  constant  field  in  die  iron  while  preventing 
saturation  in  the  highest  field  areas.  Details  of  the  finite  element  modeling  are  covered  in 
a  later  secticm  of  this  report 

Ferrous  Conductors 

In  this  machine,  iron  bars  are  used  to  carry  the  radial  magnetic  flux  as  well  as  the 
axial  load  current  in  both  the  inner  and  outer  rotor  active  regions.  The  advantage  of  this 
approach  is  diat  the  entire  volume  of  the  active  region  is  available  for  both  flux  and 
current  transport,  since  the  conductors  are  ferrous.  In  a  slotted  design,  the  functions  of 
current  and  flux  transport  are  separate  with  approximately  one  half  of  the  active  region 
volunre  available  for  each  function.  Because  of  this  combination  of  functions,  the  active 
region  of  the  iron  bar  design  is  smaller,  thereby  reducing  the  overall  machine  size.  The 
higher  resistiviQr  of  iron  is  mostly  compensated  for  by  the  increased  cross  sectional  area 
available  for  conductor  and  by  the  fact  that  the  active  region  is  ^pioximately  one  half  the 
length  of  the  active  region  in  a  slotted  design. 

Flux-Load  Current  Interactions 

There  are  interactions  between  the  load  current  and  radial  magnetic  inthiction  that 
must  be  analyzed  and  accounted  for  in  this  design.  The  axial  load  current  in  the  bars 
generates  a  circumferential  as  well  as  an  internal  magnetic  field  intensity  that  effectively 
reduces  the  pooieability  of  the  iron  bars.  This  occurs  because  the  magnetic  field 
intensity  in  the  material  is  the  vectOT  sum  of  the  radial  field  intensi^  from  the  coils  and 
the  internal  and  circumferential  field  intensities  generated  by  the  load  current  The  bars 
are  separated  fiom  each  other  by  G-10  (glass-epoxy  laminate)  insulators.  These  provide 
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a  non  ferrous  break  which  reduces  the  circumferential  field  produced  by  the  load  cunent 
This  combination  of  fields  proved  to  be  difficult  to  model  in  two  dimensions.  Due  to  the 
geometry  of  the  active  region  and  direction  of  the  fields,  unique  boundary  conditions  are 
required  in  the  finite  element  model.  The  modeling  was  finally  done  using  an  in-house 
developed  finite  element  program,  that  could  utilize  the  required  boundary  conditions, 
interactively  with  a  commercial  finite  element  package.  This  approach  allowed  us  to 
accurately  model  the  field  interactions  in  the  iron  bars.  This  is  described  in  detail  in  Ref. 
5 

Shielding 

In  a  ncmnally  conducting  design  such  as  this,  the  shielding  serves  several 
functions.  Rrst  is  the  minimization  of  die  external  stray  field.  Die  second  is  providing  a 
low  reluctance  path  for  the  main  wcaking  flux.  The  low  reluctance  path  is  required 
because  of  the  limited  number  of  ampere  turns  available  from  normally  conducting  field 
coils.  Providing  a  large  enough  coil  to  drive  the  flux  through  a  large  air  gap  would  have 
significant  negative  impacts  on  both  machine  size  and  efficiency.  The  flux  shield  also 
serves  as  the  motor's  outer  housing. 

Finite  Element  Models 

Finite  element  models  were  used  to  determine  the  overall  iron  dimensions  and  to 
verify  the  final  magnetic  design.  PE2D,  a  two  dimensional  nonlinear  finite  element 
package  from  Vector  Fields  was  used  for  the  final  analysis.  A  three  dimensional  finite 
element  package  was  available  in-house,  but  the  added  time  and  expense  of  deriving  three 
dimensional  models  were  beyond  the  scope  of  this  project  The  basic  iron  dimensions 
were  determined  by  running  many  models  and  iterating  the  dimensions  so  that  the  iron 
would  be  equally  loaded  throughout  the  magnetic  circuit  Once  the  overall  dimensions 
were  determined,  a  more  detailed  final  model  was  derived  to  verify  that  0.1  Weber  of  flux 
would  cross  the  active  region  of  the  bars.  The  dimensions  used  in  the  final  magnetic 
model  are  shown  in  Fig.  5. 

Although  the  geometry  of  the  motor  is  basically  axisymmetric,  there  are  some 
localized  three  dimensional  effects  such  as  the  flux-load  current  interaction  in  the  iron 
bars  and  voids  in  the  iron  due  to  construction  necessities  such  as  screw  holes  and  service 
holes.  Since  the  finite  element  package  was  only  two  dimensional,  these  three 
dimensional  effects  were  accounted  for  by  modifying  the  iron  BH  curves  at  the  places 
were  these  irregularities  occurred.  Modifying  the  BH  curves  to  account  for  the  flux-load 
interaction  was  a  complex  process  and  is  described  in  detail  in  Ref.  5.  The  BH  curves  for 
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the  voids  in  the  iron  were  more  simply  derived  by  estimating  the  percentage  of  void  in 
the  iron  material  for  qiecific  regions  and  then  modifying  the  BH  curve  for  solid  1010 
steel  fOT  those  regions.  The  formula  for  deriving  a  BH  point  for  iron  is  given  by: 

B'  =  Qis+tio(^-s))H  (1) 

were 

B  =  B  point  for  solid  1010  steel 
//  =  H  point  for  solid  1010  steel 
B'  =  Modified  B  point  with  void  in  material 
s  =  Iron  fraction  of  marerial 
^0  =  1-2566  X 10^  (B  =  Tesla.  H  =  At  /  m) 

Po  =  1.0  (B  =  Gauss,  H  =  Oersted) 

Noting  that  \l  =  B/  H  equation  (Eq.l)  can  be  rewritten  as 

B'  =  (Bs  +  H(l-j)Po)  (2) 

Six  regions  were  identified  as  having  a  lower  permeability  due  to  voids  in  the 
steel,  and  are  shown  in  Fig.  6. 

They  are: 

1:  (Mare  4, 96%)  The  area  of  the  irmer  rotor  core  just  inside  the  field  coil  which  has  slots 
for  the  field  coil  leads. 

2:  (Mate  8, 95%)  The  screw  holes  above  the  outer  rotor  bars  and  inside  the  irmer  rotor 
bars  that  attach  the  rotor  bars  to  the  steel  rotor  cylinders. 

3:  (Mate  7, 75%)  The  area  around  the  screw  heads  of  2. 

4:  (Mate  10, 65.4%)  The  area  above  the  screw  heads  in  the  outer  rotor  cylinder. 

5:  (Mate  6, 55%)  The  outer  rotor  flux  return  ring  which  has  cutouts  for  the  coppo*  of  the 
outer  rotor  to  cany  the  load  current  to  the  outer  rotor  collectors 
6:  (Mare  9, 27.3%)  The  area  below  the  screw  heads  in  the  irmer  rotor  cylinder. 

The  BH  curves  used  for  these  regions  are  tabulated  in  table  1,  and  a  flux  plot  fix)m  the 
model  solution  is  shown  in  Fig.  7.  The  three  dimensional  effects  did  not  reduce  the  total 
flux  as  much  as  was  expected.  The  total  reduction  was  only  about  5%  from  the  runs 
assuming  solid  steel.  The  final  finire  element  solution  predicted  a  total  flux  of  0.104 
Weber. 
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Fig.  6.  Penneability  Regions 


Table  1.  BH  Curves  for  1010  Steel  With  Air  Voids 


Field 

In 

Magnetic  induction  in  Gauss 

(100%  is  all  lion  and  no  Air) 

Oersted 

MATE^ 

MATE^ 

MATE=8 

MATE=7 

MATE^IO 

MATE=6 

MAT&=9 

100% 

96% 

95% 

75% 

65.4% 

55% 

273% 

0.00 

0 

0 

0 

0 

0 

2.24 

8944 

8587 

8497 

6709 

5850 

4920 

2443 

3.70 

12000 

11520 

11400 

9001 

7849 

6602 

3279 

6.30 

14000 

13440 

13300 

9158 

7703 

3827 

10.01 

15000 

14400 

14251 

11253 

9813 

8255 

22^6 

16000 

15361 

15201 

12006 

10472 

8810 

4384 

35.97 

15841 

15677 

12384 

10803 

9091 

4531 

55.06 

17000 

16322 

■■SSI 

12764 

11137 

9375 

4681 

75.95 

17500 

16803 

16629 

13144 

11471 

9659 

4833 

102.06 

18000 

17284 

17105 

13526 

11807 

9946 

4988 

133.01 

17582 

13908 

12145 

10235 

5147 

171.00 

19000 

18247 

18059 

14293 

12485 

10527 

5311 

216.45 

19500 

18729 

18536 

14679 

12828 

5481 

266.00 

19211 

19013 

15067 

13172 

11120 

5653 

336.20 

19693 

19492 

15459 

13523 

11426 

5841 

424.20 

1  21000 

20177 

19971 

15856 

13881 

11741 

6041 

480.25 

20419 

20212 

14064 

11904 

6150 

550.40 

21500 

20662 

20453 

16263 

14251 

6270 

659.02 

21750 

20906 

20695 

16477 

14453 

12259 

6417 

829.40 

■EISl 

20941 

16707 

12473 

6609 

1250.00 

22500 

21438 

17188 

15148 

12938 

7051 

1774.50 

23069 

22217 

17745 

15701 

13486 

7588 

23996 

23143 

22930 

18664 

16616 

14398 

8489 

5341.20 

26706 

25851 

25638 

21365 

19314 

17092 

11174 

10690.00 

32070 

31215 

31001 

26725 

24673 

22449 

16527 

26741.00 

48134 

47278 

47064 

42786 

40732 

38507 

32581 

42790.00 

64186 

63330 

63116 

58837 

56783 

54558 

48631 
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Fig.  7.  Magnetic  Circuit  with  Superimposed  Flux  Plot 


FIELD  CX)ILS 

Each  coil  consists  of  343  mms  of  hoUow  copper  conductor  wound  in  the  form  of 
a  solenoid,  producing  7650  ampere  turns.  Tlie  total  effective  magnetic  flux  produced  is 
0. 104  Webers.  The  coils  have  a  calculated  DC  resistance  of 0.0267  Q  and  a  calculated 
free  air  inductance  of  253  ^h.  At  the  operating  current  of  222  amps,  each  coil  will  have  a 
terminal  voltage  of  5.93  volts  and  dissipate  1.3  kW.  Coolant  flows  through  the  hollow 
conductor  to  remove  the  generated  heat  The  coils  are  connected  electrically  in  series  to 
minimize  the  current  carried  by  the  slip  rings.  The  cooling  circuits  are  in  parallel  to 
minimize  the  required  coolant  supply  pressure.  Hie  coils  were  manufactured  by  Mohler 
Technologies  Inc.  of  Boonville,  Indiana.  The  detail  design  of  the  coils  can  be  found  on 
drawing  A27- 19992-13  in  appendix  E. 

Conductor 

The  conductor  used  for  the  coils  is  square  copper,  OFHC 102  alloy,  0.218  in.  cm  a 
side.  It  has  a  round  hole  0.174  in.  in  diameter  for  coolant  flow  and  is  custom  extruded. 
Due  to  the  small  size  of  the  coils  and  limited  number  of  turns,  one  continuous  length  of 
conductor  could  be  used  for  each  (x>il,  thereby  eliminating  joints  within  the  coils. 

Insulation 

The  individual  conductors  are  wrapped  with  half  lap  turns  of  .007  in  mica 
paper/glass  tape  impregnated  with  B-staged  epoxy  resin.  The  final  conductor  insulaticm 
thickness  is  .014  in.  on  a  side.  The  coil  ID,  OD,  and  both  ends  have  ground  wall 
insulation  consisting  of  cross  ply  glass  and  B-staged  epoxy  resin  per  MIL-P-25421,  Rev 
B.  After  winding,  the  conductor  and  ground  wall  insulation  is  oven  cured.  During 
curing,  the  B-staged  epoxies  flow  together  forming  a  solid  conqmsite.  Both  the 
conductor  insulation  and  the  ground  wall  insulation  are  class  F  witii  a  test  voltage  rating 
of  1300  volts. 

Lead  Configuration 

Two  coil  leads,  consisting  of  extended  conductor,  are  provided  180°  apart  Both 
leads  are  bent  inward  so  that  they  are  at  the  same  radius  and  can  clear  the  bore  of  the 
other  coil.  An  adapter  is  soft  soldered  to  the  ends  of  the  leads.  This  adapter  provides  a 
round  surface  for  the  terminal  clamps  to  attach  to  as  well  as  a  barbed  end  for  the  coolant 
hose  connections.  The  terminal  clamps  are  bolted  to,  but  insulated  from,  the  inner  rotor. 
This  provides  mechanical  support  for  die  coil  leads. 
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ELECTRICAL  dRCUTT 

With  a  known  flux,  the  electrical  parameters  of  the  machine  can  be  determined. 
The  terminal  voltage  of  single  turn  homopolar  motor  is  given  by: 


V  =  ^ 
60 


(3) 


Where  <I>  is  the  flux  cut  by  the  lotOT  per  revolution  and  N  is  the  differential  rotor  speed  in 
revs/min.  For  this  motor  the  flux  is  0.104  webers;  therefore  the  terminal  voltage  is: 


0.104Ar 

60 


(4) 


At  the  nominal  design  speed  of  ±  1,200  RPM  this  results  in  a  open  circuit  terminal 
voltage  of  4.16  volts. 

Inner  Rotor 

The  electrical  circuit  of  the  irmer  rotor  consists  of  32  iron  bars,  piece  8,  two 
copper  transition  rings,  pieces  10  and  1 1,  and  two  current  collector  rings,  pieces  9  and  12. 
The  iron  bars  conduct  the  current  dirough  the  active  regitm  of  die  machine.  The  bars, 
which  are  trapezoidal  in  cross  section,  are  fastened  to  the  irmer  rotor  cylinder,  piece  3. 
The  bars  are  separated  fiom  each  otho’  by  G— 10  (glass-epoxy  laminate)  insulators.  In  a 
multi— turn  motor  these  would  provide  the  electrical  isolation  between  circuits.  They  also 
provide  a  non-ferrous  break  as  previously  described.  TTie  outer  surface  of  the  inner  rotor 
cylinder  is  plasma  sprayed  with  a  0.015  in  thick  coating  consisting  of  a  mixture  of 
aluminum  oxide  and  titanium  dioxide  (AI2O3  and  1102)  for  electrical  isolation.  Bolts 
with  insulated  shanks  and  insulating  washers,  pieces  60  and  63,  are  used  to  secure  the 
bars  to  the  rimer  rotor  cylinder.  The  copper  irmer  rotor  transition  tings  are  brazed  to  the 
ends  of  the  bar  assembly.  Hiese  tings  have  pockets  machined  in  the  iimo*  face  to 
accommodate  the  cooling  tube  U-bends,  piece  57  and  have  tapped  holes  on  the  outer 
faces  to  which  the  current  collector  disk  bolt  The  mating  faces  of  die  transition  ring  and 
the  collector  disk  are  gold  plated  to  minimize  contact  resistance.  The  cross  sectional  area 
of  the  bars  is  50.8  square  inches.  This  results  in  a  current  density  of  688  amp^square 
inch  at  35,0(X)  amps  and  1967  amps/square  inch  at  100,(XX)  amps. 
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Outer  Rotor 

The  electrical  circuit  of  the  outer  rotor  consists  of  36  iron  bars,  piece  16,  two 
coppo*  current  collector  chaiuiels,  pieces  14  and  19,  the  copper  outer  rotor  outer  drum, 
piece  17,  and  two  tominal  collector  disks,  pieces  15  and  18.  The  bar  assembly  has  a 
copper  cunent  collector,  piece  19,  brazed  on  one  end  and  the  outer  rotor  outer  drum, 
piece  17,  brazed  to  the  other  end.  This  assembly  is  insulated  with  glass-epoxy  on  its  OD 
and  is  bolted  te  the  ID  of  the  outer  rotor  cylindo*,  piece  22,  with  insulated  bolts  and 
washers,  pieces  59  and  63.  The  service  end  terminal  collector  disk,  piece  15,  is  bolted  to 
the  end  of  the  outer  rotOT  service  end  collector,  piece  14,  and  the  drive  end  terminal 
coUectcH’  disk  is  bolted  to  the  end  of  the  outer  rotor  outer  drum,  piece  17.  The  mating 
faces  are  gold  plated  as  on  the  iimer  rotor. 

Terminals 

Two  terminal  collectors ,  pieces  26  and  27  are  fitted  between  the  motor  housing, 
piece  37,  and  the  service  end,  end  plate  housing,  piece  39.  The  collectors  are  insulated 
from  each  other  and  from  the  housing  with  bonded  on  G-10  (Glass-q>oxy  laminate). 

Each  collector  has  4  input  terminal  ears,  that  protrade  through  the  terminal  housing,  piece 
38,  that  are  the  motor’s  input  terminals. 

CURRENT  COLLECTORS 

The  collectors  used  in  this  machine  are  close  clearance  liquid  metal  braid  type 
collectors.  These  collectors  have  only  been  used  in  single  rotating  applications  at 
CDNSWC  in  the  past  There  are  some  fluid  flow  and  stability  questions  concerning 
operation  of  these  collectors  in  a  contrarotating  mode  that  were  addressed  as  part  of  this 
project  The  collector  development  is  being  approached  both  theoretically  and  by  test  rig 
evaluations  of  candidate  current  collectors.  The  results  of  these  evaluations  will  be  used 
to  select  the  final  current  collector  configuration  for  both  the  single  and  contrarotating 
collectors  in  this  design.  Details  of  these  studies  can  be  found  in  Ref.  6 

COOLING  SYSTEM  DESIGN  AND  ANALYSIS 
Inner  Rotor 

The  inner  rotor  bar  set  is  the  critical  portion  of  the  inner  rotor  from  a  thermal 
perspective.  Appendix  A  documents  the  detailed  thermal  analysis  performed  tm  the 
inner  rotor  bars. 
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Selecting  coolineiiole  diameter  and  number  of  parallel  paths 

The  criterion  for  selecting  the  diameter  for  the  cooling  hole  was  to  choose  the 
minimum  cooling  hole  size  that  would  adequately  cool  the  bars  while  preserving  the 
maximum  iron  cross  section  to  carry  the  current  This  also  minimizes  the  power  loss  in 
the  armature  circuit  Figs.  8-13  display  how  the  film  drop  and  fluid  changes  in 
temperature  are  affected  by  varying  the  diameter  of  the  cooling  passage.  The  figures 
show  the  results  for  40  and  80  psi  pressures  as  well  as  different  numbers  of  parallel  paths. 
Increasing  the  diameter  decreases  the  fluid  temperature  change  at  a  higher  rate  than  the 
film  drop  temperature,  because  die  film  drop  varies  linearly  with  diameter  (surface  area) 
and  the  change  in  fluid  temperature  varies  with  the  square  of  the  diameto'  (cooling  hole 
cross  sectional  area). 

The  graphs  also  shows  how  the  temperature  changes  are  affected  by  the  number 
of  parallel  paths  in  the  system.  For  one  parallel  path,  the  coolant  must  make  a  series  run 
through  all  the  bars  before  returning  to  the  heat  exchanger.  By  increasing  the  numba*  of 
cooling  paths,  the  change  in  fluid  temperature  (less  heat  to  be  removed  per  cooling  loop) 
declines  more  than  the  film  drop  (less  heat  to  be  removed  per  cooling  loop  but  surface 
area  is  also  decreased).  Fig.  14  shows  how  the  film  drop  tempoature  varies  with 
diameter  and  pressure. 

Based  on  the  above  information,  a  diameter  of  .25"  with  2  parallel  paths  was 
chosen  to  cool  the  rotor  bars.  This  diameter  size  provided  a  good  flow  rate,  used 
miiumum  iron  bar  area,  and  had  the  ptopa*  U  bend  radius  for  construction  purposes.  Fig. 
15  displays  the  film  drop  tenqrerature  verses  pressure  and  flow  for  2  parallel  paths  and 
.25"  diameter.  The  film  drop  is  the  driving  parameter  in  determirung  rotor  bar 
temperature. 

Outer  Rotor 

The  outer  rotor  bar  cooling  analysis  is  very  similar  to  the  itmo'  rotor  bar  analysis. 
The  only  relevant  difference  is  that  the  outer  rotor  has  36  bars  each  of  a  smaller  cross 
section  resulting  in  a  greater  cooling  passage  length.  This  results  in  increased  pressure 
drop  as  well  as  increased  surface  area.  However,  the  heat  generation  in  the  outer  rotor  is 
the  same  as  the  timer  rotor  since  the  total  current,  area,  and  active  length  for  both  rotors 
are  the  same.  Appendix  B  shows  the  differences  between  the  timer  and  outer  rotor 
calculations. 

Fig.  16  shows  the  film  drop  of  the  outer  bars  for  various  flow  rates  at  2  parallel 
passes  and  a  diameter  of  .25".  The  film  drop  for  the  onto*  rotor  is  a  few  degrees  lower 
than  for  the  timer  rotor  under  the  same  conditions.  Although  the  pressure  drop  is  higher 


22 


CARDIVNS  WC-TR-8 1-94/25 


to  the  greater  number  of  bars,  the  surface  area  also  increases,  and  it  has  a  larger  effect 
on  the  film  drop  than  does  the  pressure  drop.  The  change  in  bulk  fluid  temperature 
increases  for  the  outer  rotor  compared  the  inner  rotor  results,  but  it  is  still  small  in 
magnitude  when  compared  to  the  film  drop. 
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Coolant  Hole  Diameter  (in) 

Fig.  8.  AT  Vs.  Diameter,  N=  1, 40  psi 


Coolant  Hole  Diameter  (in) 


Fig.  9.  AT  Vs.  Diameter,  N=  1, 80  psi 
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Fig.  12.  AT  Vs.  Diameter,  N=  4, 40  psi 


Coolant  Hole  Diameter  (in) 

Fig.  13.  AT  Vs.  Diameter,  N=  4,  80  psi 
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Fig.  14.  Inner  Rotor  Film  Drop  V: 
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Terminals 

Each  of  die  terminal  collector  rings  are  cooled  by  four  soldered-4n  5/16  in  copper 
tubes.  Since  these  are  completely  outside  the  cover  gas  region  of  the  motor,  water  will  be 
used  as  the  coolant  Hie  tubing  provides  70  in^  of  surface  area  per  ring  which,  when 
water  cooled,  should  be  more  than  adequate  the  anticipated  heat  load.  No  detailed 
thermal  analysis  of  the  terminal  cooling  was  performed. 

Field  coils 

Appendix  C  documents  the  thermal  analysis  procedure  used  in  designing  the  field 

coils. 

Equations  la  through  7b  of  Appendix  C  were  analyzed  to  determine  which 
combination  of  parameters  resulted  in  the  lowest  change  in  temperature  for  the  cooling 
fluid.  The  change  in  fluid  tempoxiture  was  the  driving  parameter  in  determining  the 
cooling  hole  size  and  flow  rates. 

Fluids 

Figs.  17  and  18  show  the  diffnences  in  thermal  performance  (ATfl  )  between 
Coolanol™  20  and  Coolanol™  40  heat  transfer  fluids  at  20  and  60  psi  of  pressure.  Since 
the  Coolanol™  20  produces  approximately  25  *F  less  tempo'ature  change  of  the  fluid 
than  the  Coolanol™  40  (due  to  its  lower  viscosity),  it  is  the  fluid  chosen  to  cool  the  coils. 

Geometry 

Both  the  square  and  round  cooling  holes  were  evaluated  over  a  range  of  cross 
sectional  areas  for  20  and  60  psi  pressures  as  shown  in  Figs.  19  and  20.  Fot  a  fixed 
current,  the  electrical  losses  are  dependent  on  the  cross  sectional  area  of  conductor,  but 
the  cooling  capacity  is  dependent  on  die  flow  rate,  which  in  turn  is  dependent  on  the  hole 
area.  When  the  area  of  the  hole  is  small  compared  to  the  area  of  die  conductor, 
incremental  increases  of  the  hole  dimension  yield  a  large  percentage  increase  in  flow  area 
and  a  small  percentage  decrease  in  conducts  area.  As  hole  dimension  becomes 
significant  compared  to  the  overall  dimensions  of  the  conductor,  incremental  increases  of 
the  hole  dimension  result  in  small  percentage  increases  in  the  coolant  flow  area  but  large 
percentage  decreases  in  the  conductor  area.  Since  the  round  hole  offers  better  thermal 
performance  due  to  the  higher  flow  rate  (9.4%  )  it  provides  for  the  same  cross  sectional 
area  and  pressure  drop,  the  round  hole  was  chosen  for  the  geometric  shape  of  the  cooling 
hole.  The  optimum  diameter  for  the  round  hole  (based  on  lowest  temperature  change  of 
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the  fluid)  is  found  to  be  d  =  0.185"  as  shown  in  Fig.  21.  For  reasons  of  wall  thickness 
and  availability  of  tubing  a  hole  size  of  0.174"  was  chosen.  The  difference  in 
performance  is  not  significant  as  can  be  seen  in  Fig.  21.  The  analysis  also  revealed  that 
the  optimum  diameter  is  only  a  function  of  geometric  parameters,  and  is  not  dependent  on 
fluid  properties  or  the  magnitude  of  the  heat  load.  A  diameter  smaller  than  optimum 
results  in  a  higher  temperature  rise  due  to  the  low  flow  rate  for  a  given  pressure  drop. 
Conversely,  a  larger  than  optimum  hole  will  produce  a  higher  temperature  rise  of  the 
fluid  since  the  reduction  in  copper  area  increases  die  I^R  losses  more  rapidly  than  the 
increase  in  coolant  flow  rate  can  compensate  for. 

The  change  in  fluid  tempoature  (ATfi )  is  dependent  on  the  hole  area,  while  the 
film  drop  (ATfd)  is  dependent  on  the  surface  area  of  the  cooling  hole.  Fig.  22  ^ows  the 
film  drop  temperature  difference  (round  hole)  versus  diameter  for  various  pressures.  For 
a  diameter  of  d  =  .174"  the  largest  film  drop  is  8.5  *F  at  20  psi  (lowest  pressure 
considered).  Since  the  film  drop  is  small  in  magnitude  compared  to  ATfl  thoe  is  no  need 
to  be  concerned  about  this  parameter. 

Based  on  equations  8  through  12  of  Appendix  C  Figs.  23-25  show  the 
relationship  between  pressure  and  flow  rate  and  how  they  affect  the  change  in 
ten:q)erature  of  the  fluid. 
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Fluid  A  T  (°F)  Fluid  AT  (°F) 


Fig.  17.  Fluid  AT  Vs.  Hole  Dia.  for  Coolanol  20  and  40  at  20  psi 


Fig.  18.  Fluid  A  T  Vs.  Hole  Dia.  for  Coolanol  20  and  40  at  60  psi 
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Fig.  19.  Ruid  AT  Vs.  Coolant  Hole  Area  at  20  psi 
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Fig.  20.  Fluid  AT  Vs.  Coolant  Hole  Area  at  60  psi 
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Rg.  21.  Ruid  AT  Vs.  Coolant  Hole  Diameter 
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Coolant  Hole  Diameter  (in) 


Fig.  22.  Film  Drop  AT  Vs.  Diameter  at  20, 40, 60,  and  80  psi 


Fig.  23.  Huid  AT  Vs.  Flow  Rate 
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STRUCTURAL  DESIGN 

Stress  and  deflection  analyses  were  performed  on  the  rotating  elements  of  the 
motor.  This  included:  (1)  manual  calculations  using  traditional  "textbook"  methods  to 
provide  an  adequate  design  starting  point  and  (2)  a  finite  element  analysis  during  the 
detail  design  stage  when  the  final  geometry  was  more  firmly  established.  The  following 
surrunarizes  the  assumptions  and  results  of  the  first  analysis  while  appendix  D  presents 
the  results  of  the  second  analysis  which  determined  the  final  detail  design. 

While  the  motOT  has  a  nominal  rating  of  175  hp  at  ±  1200  rpm  (87.5  hp  per  shaft), 
the  overload  rating  of  6(X)  hp  was  used  (300  hp  per  shaft)  for  the  analysis.  As  the  method 
for  loading  the  motor  was  unknown  at  the  time  of  the  analysis,  a  Browtung  "Poly- V" 
multi-rib  V-belt  drive  was  selected  as  a  worst  case  (largest  overhimg  load)  method  for 
transfer  of  power  to  the  load  device.  Combiiung  sheave  weight  and  belt  tension  forces, 
an  overhung  load  of  400  lb  with  an  assumed  moment  arm  of  24  inches  was  used  for  the 
analysis.  Rotor  weights  were  conservatively  estimated  based  on  the  weight  of  solid 
copper  resulting  in  timer  and  outer  rotor  weights  of  526  lb  and  739  lb,  respectively.  The 
analysis  was  conducted  using  the  smallest  diameto’  portion  (2.953  inch  O.D. )  of  the 
timer  shaft  as  the  structural  member.  All  analyses  were  performed  using  304  stainless 
steel  shafting  with  a  yield  strength  of  35,000  psi  and  a  fully  modified  endurance  limit  of 
18,100  psi. 

The  stress  analysis  was  performed  using  three  methods.  The  first  method  was  for 
static  loading  conditions  and  resulted  in  a  maximum  shear  stress  of  3650  psi  and  a  factor 
of  safety  of  4.8.  The  second  method  was  fOT  revised  bending  and  steady  torsion  and 
resulted  in  an  alternating  stress  of  3790  psi  and  a  factor  of  safety  of  4.8.  The  third 
method  was  the  Soderberg  approach,  which  is  very  conservative  and  appears  to  be  the 
basis  for  the  Navy  propulsion  shafting  design  data  sheet  The  Soderberg  approach 
resulted  in  a  required  minimum  shaft  diameter  of  1.91  inches  and  a  factor  of  safety  of  3.6. 
Applying  a  stress  concentration  factOT  of  U  for  the  shaft  step  reduced  the  three  factors  of 
safety  to  3.2, 3.2,  and  2.4,  respectively. 

The  deflection  analysis  predicted  radial  deflections  of 0.0021  and  0.0026  inches 
for  the  timer  rotor  service  end  and  drive  end,  respectively.  Using  the  ratio  of  timer  to 
outer  rotor  loads  and  moments  of  inotia,  an  outer  rotor  maximum  deflection  of  less  than 
0.0006  inches  is  predicted.  Thus,  an  timer  rotor  to  outer  rotor  maximum  net  deflection  of 
approx.  0.002  inches  is  anticipated.  The  predicted  outer  rotor  to  motor  casing  maximum 
net  deflection  is  0.0006  inches. 
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These  initial  stress  and  deflection  analyses  utilize  conservative  assun:q)tions  in 
regards  to  loading  and  structural  stiffness.  The  outer  shaft  is  much  stronger  than  the  inner 
shaft  and  the  center  portion  of  the  inner  and  outer  rotor  are  much  stronger  than  the  necked 
down  shaft  sections.  Therefme,  general  stresses  are  not  considered  to  be  a  problem. 
Similarly,  while  maximum  net  radial  deflections  of 0.002  inches  were  calculated,  the 
actual  deflection  will  probably  be  less  than  0.001  inches. 

BEARINGS  AND  GAS  SEALS 
Bearings 

Due  to  the  necessity  of  placing  bearings  between  the  inner  and  outer  rotors  and 
the  desire  to  keep  the  machine  as  compact  as  possible,  bearings  of  the  smallest  cross 
section  that  had  an  adequate  load  rating  were  selected.  To  prevent  eddy  currents  through 
the  bearings,  the  surface  of  die  shafts  are  coated  with  a  plasma  deposited  ceramic  undo* 
the  inner  race  of  the  bearings. 

Drive  end 

The  bearings  selected  for  the  drive  end  are  deep  grove  ball  bearings  which  are 
grease  prelubricated  and  shielded,  pieces  49  and  SO.  Tbese  bearings  provide  both  axial 
and  radial  support  and  location  for  the  rotors.  The  pairs  of  bearings  will  be  shimmed  to 
minimize  play. 

Service  end 

For  the  service  end,  double  row,  cylindrical  roller  bearings  were  chosen  to  provide 
radial  support  as  well  as  to  accommodate  axial  thermal  growth,  and  are  of  JIS  Oass  4 
precision,  grease  prelubricated  and  shielded,  pieces  45  and  46. 

Seal  Assembly 

Gas  seal 

Due  to  the  reactive  nature  of  NaK  a  dry,  oxygen  free  cover  gas  must  be 
maintained  in  the  machine  at  a  pressure  several  PSI  above  atmosphere.  This  necessitates 
the  use  of  high  performance  gas  seals  on  both  shafts,  pieces  43, 44, 47,  and  48. 

Mechanical  face  seals  provide  both  the  required  performance  and  low  leakage  rate.  A 
carbon  face  running  on  a  AI2O3  coating  on  the  shoulder  has  worked  well  in  similar 
applications  in  the  past  and  will  be  utilized  here.  The  AL2O3  coating  provides  an 
insulating  break  to  prevent  eddy  currents  through  the  seals.  Most  carbons  depend  on 
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atmospheric  moisture  to  form  a  lubricating  film.  Since  these  seals  operate  in  a  very  dry 
environment  a  special  grade  of  carbon  (Pure  Carbon  P5N)  is  used  to  prevent  rapid  wear. 

Exclusion  seal 

To  prevent  contamination  of  the  face  seals  with  bearing  lubricant,  exclusion  seals, 
pieces  103, 104,  and  105,  are  located  between  the  face  seals  and  the  bearings.  These 
seals  are  of  a  simple  lip  design.  They  utilize  a  self  lubricating  plastic  sealing  element,  in 
a  metal  housing,  miming  on  hardened  steel  sleeves,  pieces  71, 72, 77,  and  79,  <hi  the 
shafts. 

SERVICE  MODULE 
SUP  RINGS 

A  slip  ring  assembly,  to  supply  current  to  the  field  coils,  is  attached  to  the  service 
module.  This  is  of  conventional  design  consisting  of  copper  slip  rings  and  silver  graphite 
bmshes.  The  copper  terminal  rods  of  the  inner  rotor  extend  through  the  slip  ring 
assembly  so  that  connections  to  the  slip  rings  can  be  made.  The  slip  ring  assembly  is  a 
commercial  unit  produced  by  Fabricast  Inc.,  So.  El  Monte,  CA  91733 

COOLANT  SEALS 

The  service  module  contains  a  combination  of  mechanical  face  seals  and  Up  seals 
to  provide  coolant  flow  paths  to  and  from  the  inner  and  outer  rotor.  The  face  seals 
provide  both  the  required  performance  and  low  leakage  rate.  A  carbon  face  running  on  a 
AI2O3  coating  on  the  shoulder  is  utiUzed.  The  Up  seals  utilize  a  self  lubricating  plastic 
sealing  element,  in  a  metal  housing,  running  on  a  hardened  steel  sleeve  on  die  shaft  The 
Up  seals  are  used  in  the  less  critical  locations  in  the  service  module  where  some  leakage 
can  be  tolerated. 

ROTARY  UNION 

A  rotary  union  is  fitted  on  the  end  of  the  service  module,  outboard  of  the  sUp  ring 
assembly,  to  supply  coolant  to  the  machine. 

INSTRUMENTATION 

Sensors 

Thermocouples  are  fitted  to  measure  the  temperature  of  coolant,  iimer  and  outer 
rotor  bars,  inner  and  outer  rotor  bearings,  and  current  coUectors. 
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Telemetry 

Provisions  were  made  in  the  design  to  allow  both  inner  and  outer  rotor  to  be  fitted 
with  telemetiy  systems  in  the  future.  They  will  consist  of  a  multiplexing  unit  and  a 
transmitter  on  each  rotm*  which  will  transmit  the  theimocouple  data  to  a  fixed  receiver. 

COVER  GAS  SYSTEM 

The  space  between  the  inner  and  outer  rotor  and  between  the  outer  rotor  and 
housing  is  supplied  with  high  purity  nitrogen  cover  gas  to  prevent  reaction  of  the  sodium- 
potassium  used  in  the  cuiient  collectors  with  oxygen  tx  wato*  vtqxir.  A  pre-purified 
grade  of  nitrogen  is  generally  supplied  from  K  cylinders  and  is  regulated  to  several  psig. 
The  cover  gas  is  supplied  to  both  ends  of  the  motor  to  minimize  flow  through  the  motor 
which  would  increase  the  reaction  with  the  liquid  metal. 

ASSEMBLY  PROCEDURES 

Assembly  will  begin  with  the  inner  rotor  which  will  be  fitted  with  the  bar  set  and 
the  field  coils.  Coolant  connections  will  be  made  with  flexible  hose  to  the  passages  in  the 
shaft  Coil  electrical  connections  will  be  made  with  split  clanq)s  and  rigid  copper  rods 
inteicotmecting  with  the  terminal  rods  in  the  shaft  The  outer  rotor  bar  set  is  semi¬ 
permanently  mounted  intonal  to  the  outer  rotOT  cylinder.  Both  inner  and  outer  rotor 
current  collector  rings  and  disks  will  be  bolted  in  place  sequentially  as  the  outer  rotor 
assembly  in  built  up  around  the  itmer  rotra*.  The  flux  return  rings  will  dien  be  bolted  in 
place.  Bearings  and  seals  will  be  installed  in  the  outer  rotor  stub  shafts  which  will  tiien 
be  bolted  to  tiie  outer  rotor  assembly.  Bearings  and  seals  will  be  fitted  in  the  outer  rotor 
bearing  carriers  which  will  then  be  installed  over  the  outa*  rotor  shafts.  These  bearings 
and  seals  need  not  be  disturbed  to  accomplish  the  level  of  disassembly  required  for 
cleaning  or  other  minor  maintenance.  The  assembly  will  be  inserted  in  the  housing  and 
the  various  termiiud  collectOT  components  installed.  The  coolant  connections  for  the 
outer  rotor  components  will  be  made  with  flexible  hose  to  the  passages  in  tiie  stub  shafts. 
Installation  of  the  two  outer  rotor  bearing  housings  completes  the  basic  machine 
assembly. 


PREDICTED  MACHINE  PERFORMANCE 
Output  power  calculations  for  a  homopolar  motor  are  rather  straight  forward.  The 
open  circuit  terminal  voltage  is  only  a  function  of  machine  flux  and  rotor  speed.  The 
internally  generated  torque  is  a  function  of  motor  flux  and  load  current  The  actual 


39 


CARDIVNSWC-TR-^1-94/25 


terminal  voltage  and  output  torque,  however,  are  affected  by  the  motor’s  internal  losses. 
These  losses  can  be  divided  into  three  basic  categories. 

Resistive  voltage  drops  in  the  intonal  conductors,  joints,  and  current  collectors 
add  to  the  generated  back  EMF  thereby  increasing  the  terminal  voltage  of  the  motor.  The 
increase  in  terminal  voltage  is  linear  with  load  current  The  power  loss  due  to  the  motor's 
internal  resistance  is  proportional  to  the  square  of  die  load  current  Both  are  primarily 
independent  of  motor  speed  and  flux. 

Frictitxial  drag  subtracts  from  the  actual  internally  genoated  torque  and  reduces 
the  output  torque  produced  by  the  machine.  IrictitHial  drag  from  the  tolling  element 
bearings  and  seals  is  relatively  ccmstant  and  does  not  change  with  motor  speed  or  load. 
Windage  drag  is  a  function  of  motor  speed  and  geometry  and  is  ccxnpletely  independent 
of  load. 

The  losses  in  the  liquid  metal  current  collectors  consist  of  resistive,  frictional, 
fluid  shear,  eddy  current,  and  magneto  hydrodynamic  type  losses.  These  losses  are  a 
complicated  function  of  collector  speed,  load  current,  geometry,  magnetic  field,  and  the 
fluid  properties  of  die  liquid  metaL 

Tables  3  and  4  detail  the  predicted  performance  of  die  motor  at  35,000  and 
100,000  ampraes  load  current  respectively. 

Table  2.  Predicted  Motor  Performance  at  35,000  ampoes 


185  hp 

Shaft  speed 

±1,200  RPM 

4051b-ft 

Effective  magnetic  flux 

.104  Webers 

Back  EMF  (at  ±  1,200  RPM) 

4.16  volts 

Terminal  voltage  (at  35,(X)0  amps) 

4.28  volts 

Input  current 

35,(X)0  amps 

Estimated  efficiency 

92% 
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Table  3.  Predicted  Motor  Performance  at  100,000  amperes 


545  hp 

Shaft  speed 

±1,200  RPM 

1,192  Ib-ft 

Effective  magnetic  flux 

.104  Webers 

Back  EMF  (at  ±  1,200  RPM) 

4.16  volts 

Terminal  voltage  (at  1(X),(XX)  amps) 

4.50  volts 

Input  current 

1(X),0(X}  amps 

Estimated  efficiency 

90% 

CONCLUSIONS 

From  a  basic  point  of  view  a  contrarotating  hcnnopolar  motOT  is  simply  one  in 
which  the  stator  is  allowed  to  rotate  and  produce  powCT.  From  the  outset  we  realized  that 
this  would  be  a  difficult  design  challenge.  It  turned  out  to  be  much  more  difficult  than 
we  anticipated.  Providing  the  required  magnetic  flux,  adequate  cross-section  for  the 
current,  sufficient  coolant  flow,  and  providing  the  structure  required  to  transmit  the 
torque  did  not  turn  out  to  be  the  most  difficult  part  of  the  design.  The  real  complexity 
turned  out  to  be  the  multitude  of  details.  All  the  rotating  members  must  be  supplied  with 
coolant  This  required  a  number  of  connections,  passages,  seals,  and  tenq>erature  sensors 
for  both  the  inner  and  outer  rotOT.  Both  rotors  must  be  fitted  with  bearings  and  seals. 

The  sizes  of  suitable,  available  bearings  had  a  significant  impact  on  shaft  sizes.  A  great 
deal  of  effort  also  went  into  the  assembly  /  disassembly  procedures  since  this  machine 
does  not  have  a  split  outer  rotor  or  outer  housing.  Although  simple  in  concept,  the 
contrarotating  homopolar  design  is  not  simple  to  execute.  Although  it  took  considerable 
time  in  the  design  phase,  we  were  able  to  come  up  with  a  design  that  provided  for  the 
many  unique  requirements  of  diis  motor.  At  this  time,  motor  construction  is  complete 
and  assembly  is  underway.  Based  on  our  experience  in  ccHistruction  and  assembly  of  this 
motor  the  basic  design  appears  viable  and,  with  some  changes,  suitable  for  largo*  motors. 
In  fact,  a  larger  motor  presents  fewer  difficulties  because  it  removes  the  volume 
constraints  of  a  small  motor. 
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APPENDIX  A 

INNER  ROTOR  THERMAL  DESIGN 


The  inner  rotor  bar  set  is  the  critical  portion  of  the  inner  rotor  from  a  thermal 
perspective.  This  section  documents  the  thermal  analysis  performed  on  the  inner  rotor 
bars.  The  thermal  analysis  is  presented  in  three  sections  that  include  design  constraints, 
assumptions,  and  the  results. 

DESIGN  CONSTRAINTS 

The  geometry  of  the  inner  rotor’s  bars  was  determined  by  magnetic,  electrical, 
and  space  considerations.  The  length  of  the  bar  was  driven  mainly  by  magnetics.  Height 
of  the  bar  was  driven  by  space  in  the  machine.  Width  and  number  of  bars  were 
determined  by  electrical,  cooling,  and  insulation  charac^stics.  There  are  32  bars  on  die 
mner  rotor  that  each  have  a  length  of  5  inches,  and  liquid  coolant  is  passed  axially 
through  the  center  of  each  ircm  bar. 

ASSUMPTIONS 

Based  on  the  thermal  analysis  accomplished  (Hi  die  machine's  field  coils,  the 
thermal  analysis  performed  on  the  rotor  vdll  also  use  a  round  (xioling  hole,  aatunift 
constant  properties  over  the  machine's  temperature  range,  and  use  Coolanol™  20  as  the 
cooling  flui(L  Ctxilanol™  20  displays  g(xxi  thermal  characteristics,  has  low  electrical 
conductivity,  and  does  not  react  with  liquid  metal  (NaK).  Since  this  motor  rotates  at  low 
speeds  and  has  a  small  rotcw  radius,  it  was  assumed  that  the  centrifugal  forces  acting  on 
the  (xxiling  fluid  would  have  no  significant  effect  on  its  performance. 

RESULTS 

The  results  consist  of  5  sections  of  calculations  consisting  of  evaluating  geometry, 
heat  generation  of  the  bars,  head  loss  and  flow  rates,  (xinvection  (xieflicient,  and  film 
drop  and  fluid  change  in  temperature.  The  overall  objet^ve  was  to  determine  the 
differences  between  the  input  and  ouq>ut  temperatures  of  the  fluid  (ATfl)  and  the 
differences  between  the  average  temperatures  of  the  (X)pper  and  fluid  (ATfd  -  film  drop), 
verses  the  required  flow  rate  and  pressure  drop. 


LIST  OF  ABBREVIATIONS 
Acool  =  total  (xwling  area  -  (in^) 

Ag  =  surface  area  of  ccmling  hole  ft^ 
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Cp  =  specific  heat  Btu/llm  -*F 
d  =  dimension  of  the  cooling  hole  (in) 

f  =  friction  factor  across  the  coils  inner  surface  (dimensionless  number) 

gc  =  gravitational  constant  32.2  Ibm  fr/lbf  sec^ 

h  =  convection  coefficient  Btu/hr-ft-*F 

k  =  thermal  conductivity  Btu/hr.- ft  -T 

L  =  length  of  the  conducts  ft 

Le  =  equivalent  length  ft 

m  =  mass  flow  rale  Ibmfrnin 

N  =  numbo*  of  parallel  paths 

Nu  =  Nusset  numbo*  (dimensionless  number) 

q  =  flow  rate  gpm 

Q  =  heat  generation  Btu/min 

P  =  wetted  perimeter  in 

Pr  =  Prandd  number  (dimensionless  number) 

r  =  resistivity  of  copper  ohm-in 

Re=  Reynolds  number  (dimensionless  number) 

s  =  height  and  width  of  conductor  -  in 

a  thermal  diffiisivity  ft^/hr 

AP  =  pressure  drop  (psi) 

ATf(j  =  film  drop  ((ATfd  =Ts-Tni) 

ATfl=  temperature  change  of  fluid  *F(ATfl  =Touq>ut-1'input  of  fluid) 

Ts  =  surface  temperature  of  the  conducting  metal 

Tm  =  mean  terr^rature  of  the  fluid  at  a  cross  section  in  die  flow 

p  =  density  Ibm/ft^ 

|i.  =  dynamic  viscosity  Ibm/ft-sec 

TABLE  A-1.  Fluid  Properties  of  Coolanol™  20  at  lOOT 


Properties 

Units 

Coolanol™  20 

density 

Ibm/ft^ 

55.2 

iynamic  viscosity 

Ibm^t-sec 

.0011886 

Btu/lbm  -*F 

.47 

thermal  conductivity 

Btu/hr.-ft  -’F 

.067 
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ROUND  HOLE  GEOMETRY 

Inner  rotor  -  32  bars 
d  =  dimension  of  the  cooling  hole  (in) 

L  =  actual  length  of  the  bars  in  a  parallel  path  -  5/12*  *  32/N  =  13,33/N  (ft) 
(bar  length  =  5  inches) 

32/N  =  number  of  bars  /  parallel  path 

32/N  - 1  =number of U- bends /per parallelpath 

Equivalent  length  of  head  loss  to  that  of  a  straight  pipe  for  U  -  bends 

Lc  (180®  U  -  bends)  =  50d/12 

Equivalent  length  of  head  loss  to  that  of  a  straight  pipe  -  total 
Lc  (U  -  bends  +  bars)  =  32/N(5/12)  +  50(1/12(32/N  - 1) 

Cross  sectional  area  of  the  cooling  hole 

Acool  =  cooling  cross  sectional  area  -  Jcd2/4(144)  (1) 

=  .00545d2ft2  (1.1) 

Surface  area  of  the  cooling  hole  (lengdi  times  perimeter) 

As  =  7C  dL/12N  =  Jid/12  (32  (5/12))/N 

=  3.49d/N(ft2)  (1.2) 

HEAD  LOSS  AND  FLOW  RATES 

Flow  was  assumed  to  be  fully  turbulent,  smooth  tubing,  and  incompressible  fluid, 
f  =  friction  factor  across  the  coils  inno*  surface 

(a)  Rc  =  pVd4i 

Re=  Ratio  of  die  inertia  and  viscous  forces  (Reynolds  number) 

(b)  f  =  .32/Re  *25  =  .3^-25/(V-25  (i-25p.25)  fta-  Re<  10  5  -  Blasius  correlation 

AP  =  pressure  drop 

(c)  AP  =  f  (Le/d)  (pV?/2gc)  (straight  pipe  assunption) 
substituting  (b)  into  (c) 

AP  =  .ItiLcV  1.75p.25p.75/(d  ^-^^gg)  (psi)  (2) 

substituting  for  the  fluid's  properties  into  (2) 

AP  =  .00287LcV  1-75/d  1-25  (psi)  (2.1) 

therefore  the  velocity  is 

V  =  (349APdl-25/Le>57143(f,/s)  (2.2) 

where: 

Lc  =  ft ,  V  =  ft/sec ,  d  =  in ,  p  =  Ibrnm^ ,  AP  =  psi ,  p  =  .001 1886  ft2 /sec , 
gc  =  32.2  Ibm  fl/lbf  sec^ 

(d)  m  =  pVA  mass  flow  rate 
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m=1084d2v  Ibm/hr 

(3) 

flow  rate  =  q  =  VA 

Flow  rate  per  parallel  path  on  the  iimer  rotes'  section 
q  =  2.45d2v  (gpm/parallel  path) 

Total  flow  rate  over  the  inner  rotor 

Nq  =  total  gpm 

(4) 

CONVECTION  COEFFICIENT 

Nusset  number,  and  colbum  equation 

(a) 

Nu  =  .0223  Re-8  Pr  -333  .  colbum  equatiem  where  assumptions 
(.7<Pr<160,L/D>10) 

(b) 

Pr  =  v/a  (Prandtl  number) 

Pr  =  Ratio  of  momentum  and  thermal  diffusivities 

(c) 

a  =  K/pCp  (thermal  diffusivity  (Re  >  2300) 

(d) 

h  =  NuK/d  (convection  coefficient) 

substituting  (a),  (b),  and  (c)  into  (d) 

h  =  .0223p-8V-8Cp-333K.667/(ji.467d^)  (5) 

substituting  for  the  fluid's  properties  into  (5) 

h  =  42.57  V-8  /d-2  (Btu/hr.-  ft2  -’F)  (5.1) 

where: 

V  =  ft/sec ,  d  =  in ,  p  =  lbm^t3 ,  p  =  .0011886  ft2  /sec ,  Cp  =  Btu/lbm  -T , 

K  =  Btu/hr.- ft -T,  h  =  Btu/hr.- ft2 -T 

FILM  DROP  AND  FLUID  CHANGE  IN  TEMPERATURE 
Film  drop: 

Q  =  hA,ATf(i  (ATfd  =Ts-Tin)  (mean  difference  between  the  mean  temperature 
of  the  fluid  and  the  iron  bar) 

(a)  ATfd  =  Q/AsNh 

substituting  (5.1)  &  (1.2)  &  Q  =  17,802  Btu/hr  into  (a) 

ATfd=119.8/V%8  (®F)  (6) 

whoe: 

h  =  Btu/hr.- ft2 -T ,  d  =  in ,  Q  =  Btu/min ,  ATfd  =  T 


Fluid  change  in  temperature: 

Q  =  mCpATfi  (ATfl  =  T  output  *  Tinput  of  fluid) 
(b)  ATfl=Q/(mCp) 
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substituting  (3)  &  Cp  =  .47  Btu/lbm  -*F  &  Q  =  17,802  Btu/hr  into  (b) 

ATfl  =  34.94/Nd2v  (®F)  (7) 

whoe: 

V  =  ft/sec.  Cp=  Btu/lbm -*F,  d  =  in ,  Q  =  Btu/min ,  ATfl  =  T 

The  above  equations  were  analyzed  few  determining  the  cooling  hole  diameter, 
numbCT  of  parallel  cooling  padis,  and  pressures  needed  to  adequately  cool  the  rotor  bars. 
Unlike  the  field  coil  analysis  (selecting  optimum  parameters),  the  objective  here  was  to 
chose  the  minimum  parameters  required  to  adequately  cool  the  bars.  A  minimum  number 
of  parallel  paths  is  desirable  in  ordo*  to  reduce  the  mechanical  difficulty  of  rotating  seals 
and  balancing  die  system. 
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APPENDIX  B 

OUTER  ROTOR  THERMAL  DESIGN 

The  only  changes  in  the  calculations  from  the  inner  rotor  are  noted  below: 

ROUND  HOLE  GEOMETRY 

Outer  rotor-  36  bars 
d  =  dimension  of  the  cooling  hole  (in) 

L  =  actual  length  of  the  bars  in  a  parallel  path  -  5/12*  *  36/N  =  15/N  (ft) 

3^N  =  number  of  bars  /  parallel  path 

36^  -  1  =  number  of  U  -  bends  /  per  parallel  path 

Equivalent  length  of  head  loss  to  that  of  a  straight  pipe  -  total 

Le  (U  -  bends  +  bars)  =  36/N(5/12)  +  50d/12(36/N  -  1) 

Surface  area  of  the  cooling  hole  Gengtfa  times  perimeter) 

As  =  7C  dIV12N  =  Jcd/12  (36  (5/12))/N 
=  3.93d/N(ft2) 
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APPENDIX  C 

FIELD  COIL  THERMAL  DESIGN 


The  objective  of  titis  section  is  to  document  the  thermal  analysis  procedure  used 
in  designing  the  field  coils.  The  procedure  for  the  thermal  analysis  is  broken  up  into  four 
sections  which  include,  design  constraints,  assumptions,  design  variables,  and  the  results. 

DESIGN  CONSTRAINTS 

The  coil's  geometry  was  dictated  by  magnetic,  electrical,  and  space 
considerations.  The  only  unknowns  woe  the  size  and  shape  of  the  cooling  hole  needed 
to  adequately  remove  the  heat  generated  by  the  field  cuirent  A  cross  sectional  view  for  a 
single  turn  of  the  coil  is  shown  below: 


There  are  34.5  turns  per  coil,  giving  a  total  length  of  80  ft  of  tubing  per  coil.  The 
coils  have  inner  and  outer  diameters  of  7.42"  and  9.6"  respectively  and  are  2.275"  long. 

The  coil  produces  7650  ampere  turns,  with  a  winding  cross  sectional  area  of  2.25 
in2,  thus  yielding  a  overall  current  density  of  3400  amps  per  square  inch. 

Removing  the  heat  generated  by  die  cuirent  will  be  accomplished  by  liquid 
cooling.  The  fluids  considered  were  Coolanol™  20  and  Coolanol™  40  (silicate  ester), 
because  they  exhibit  low  electrical  conductance  and  relatively  good  thermal 
characteristics  (dielectric  heat  transfer  fluid). 

ASSUMPTIONS 

During  the  thermal  analysis,  two  assumptions  were  made  in  order  to  simplify  the 
assessment.  Hie  first  assumption  established  the  fluid's  properties  anH  ct^per  resistivity 
for  the  cooling  system.  It  was  assumed  that  the  properties  (which  are  temperature 
dependent)  would  not  change  enough  to  affect  the  final  calculations.  All  die  properties 
were  taken  at  100  T  (estimated  average  temp  of  the  fluid),  and  remain  constant  through 
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the  system.  The  second  assumption  relates  to  the  pressure  drop  or  head  loss  along  the 
length  of  the  coil.  The  head  loss  was  assumed  to  be  equivalent  to  a  straight  pipe  analysis; 
that  is,  the  coil's  curvamre  was  not  taken  into  account 

DESIGN  VARIABLES 

In  the  thermal  analysis  of  the  field  coils,  there  were  two  different  variables  that 
needed  to  be  analyzed  for  the  purpose  of  determining  the  best  thermal  performance.  They 
were  the  geometry  of  the  cooling  hole  (square  or  round)  and  the  type  of  fluid  (Coolanol™ 
20  or  40).  Each  of  the  two  fluids  were  analyzed,  based  on  their  thermal  propoties,  to 
determine  which  fluid  produced  the  better  cooling  characteristics.  In  the  case  of  the 
geometry  of  the  cooling  hole,  the  square  hole  and  the  round  hole  were  compared  to  each 
other  to  determine  which  one  would  produce  die  betto'  flow  rate  fm*  the  same  area. 

RESULTS 

The  results  consist  of  6  sections  of  calculations  diat  are  comprised  of  evaluating 
the  fluids'  properties,  geometry,  heat  generation  of  the  coil,  head  loss  and  flow  rates, 
convection  coefficient,  and  film  drop  and  fluid  change  in  temperature.  The  overall 
objective  was  to  determine  the  differences  between  the  input  and  ouqiut  temperatures  of 
the  fluid  (ATfi)  and  the  differences  between  the  average  temperatures  of  the  copper  and 
fluid  (ATfd  -  film  drop),  verses  the  required  flow  rate  and  pressure  drop.  At  the 
conclusion  of  the  calculations,  the  above  assumptions  were  compared  to  a  more  realistic 
approach  to  detomine  their  validity. 

LIST  OF  ABBREVIATIONS 
Acool  =  total  cooling  area  -  (in^) 

Acu  =  conductor  cross  section  -  (in^) 

At  =  total  cross  sectional  area  -  (in^ 

As  =  surface  area  of  cooling  hole  ft^ 

Cp  =  specific  heat  Btu/lbm  -*F 
d  =  dimension  of  the  cooling  hole  (in) 

D  =  average  diameter  of  the  coil  (in) 

Dh  =  4Acool/P=bydraulic  diameter 

f  =  friction  factor  across  the  coils  inner  surface  (dimensionless  number) 
gc  =  gravitational  constant  32.2  Itnn  fl/lbf  sec^ 
h  =  convection  coefficient  Btu/hr-ft-*F 
k  =  thermal  conductivity  Btu/hr.- ft  -T 
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Jc  =  current  density  over  the  copper  amps/in^ 

Ji  =  current  density  of  the  entire  coil  amps^^ 

L  =  length  of  the  conductor  ft 

m  =  mass  flow  rate  Ibm/min 

Nu  =  Nusset  number  (dimensionless  number) 

q  =  flow  rate  gpm 

Q  =  heat  gen^tion  Btu/min 

P  =  wetted  perimeter  in 

Pf  =  packing  factor 

Pf  =  Prandd  number  (dimensionless  number) 

r  =  resistivity  of  coppo"  ohm-in 

Re  =  Reynolds  number  (dimensionless  number) 

s  =  height  and  width  of  conductor  -  in 

t  =  thickness  of  the  insulation  -  in 

vc  =  volume  of  ccmductor  in^ 

vt  =  total  volume  of  coil  in^ 

a  =  thermal  diffusivity  ft^/hr 

AP  =  pressure  drop  (psi) 

ATfd  =  film  drop  ((ATfd  =Tc-Tf) 

ATfl  =  temperature  change  of  fluid  *F(ATfl  =  T  output  *  Tinpia  of  fluid) 

p=  density  Ibm/ft^ 

p.  =  dynamic  viscosity  Ibm/ft-sec 


TABLE  C-1.  Ruid  Properties  of  Coolanol™  20  and  40  at  100®F 


Properties 

Units 

Coolanol™  20 

Coolanol™40 

densitv 

lbm/ft3 

55.2 

54.91 

Ibm/ft-sec 

.0011886 

.00402 

Btu/lbm  -*F 

.47 

.48 

thermal  conductivity 

Btu/hr.-ft  -*F 

.067 

.078 

GEOMETRY 
Round  hole: 

s  =  height  and  width  of  conductor  -  .25  (in) 
d  =  dimension  of  the  cooling  hole  (in) 
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t  =  thickness  of  the  insulation  •  .016  (in) 

At  =  total  cross  sectional  area  -  s^  (in^)  =  .0625  (in^) 

Acu  =  conductor  cross  section  -  (s-2t)2  -  jcd?/4  (in2) 

Acool  =  total  cooling  area  -  Jtd^M 

Pf  =  packing  factor  -  Acu/At  =  [(s-2t)2  -  jtd2/4]/s2  -  .7604  -  12.56<P 


Square  hole: 

At  =  total  cross  sectional  area  -  s^  (in^) 

Acu  =  wmductOT  cross  section  -  (s-2t)2  -  d^  (in^) 

Acool  =  total  cooling  area  -  d^ 

Pf  =  packing  factor  -  Acu/At  =  ([(s-2t)2  -  d2]/s2)  -  .7604  -  16d2 
Dh  =  4Acool/P=hydraulic  diameter  used  fw  correlating  turbulent  pipe  flow 
for  use  with  non  circular  geometry’s.  For  a  square  Dh  =  4d?/4d  =  d  (P  = 
wetted  poimeter). 

vc=  volume  of  conductor  (in^)  =  Acu  L  =  At  Pf  L 
L  =  length  of  the  conductor 


HEAT  GENERATION 

(a)  Q  =  J^c  r  vc  (Watts)  (heat  generation  across  the  conductor) 

(b)  Jc  =  J/Pf  current  density  over  the  conductor 
Ji  =  current  density  of  the  entire  coil 

(c)  vt  =  A^  total  volume  of  coil 

(d)  Vc  =  vtPf  volume  of  conductor 
substituting  (c)  into  (d): 

(e)  Vc  =  AtLPf 
substituting  (b)  and  (e)  into  (a): 

Q  =  (J,/Pf)2  r  (AtLPf)  =  AtUPf 
converting  to  English  units: 

(f)  Q  =  6826Jt^r  s^L/[(s-2t)2  -  jcd2/4]  Btu/min  (round  hole) 
r  =  7.87e-7  ohm-in  (resistivity  of  copper  at  100  ®F) 
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substituting  Jt  =  3400  amps^2^  s  =  .25",  t  =  .016",  and  L  =  80*  into  (f) 

Q  =  1.941/(.047524  -  .785d2)  Btu/min  (round  hole)  (la) 

Q  =  6826Jt2r  s^l7[(s-2t)2  -  d^]  Btu/min  (square  hole) 

Q  =  1.941/(.047524  -  d^)  Btu/min  (square  hole)  (lb) 

where: 

L  =  ft,s  =  in,t  =  in,Q  =  Btu/min ,  d  =  in ,  r  =  ohm-in ,  J  = 


HEAD  LOSS  AND  FLOW  RATES 

Flow  was  assumed  to  be  fiilly  turbulent,  smooth  tul^g,  and  incompressible  fluid, 
f  =  friction  factor  across  the  coils  inner  surface 


(a) 

Rc  =  pVd4i 

Re=  Ratio  of  die  inertia  and  viscous  forces  (Reynolds  number) 

(b) 

f  =  .32/Re  -25  =  .32h-25/(v.25  d.25p.25)  (RH  &  SH)  for  Re<  10  5  - 

Blasius 

correlation 

AP  =  pressure  drop 

(c) 

AP  =  f  (L/d)  (pV2/2gc)  (Horizontal  pipe  assunqition) 

substituting  (b)  into  (c) 

AP  =  1.9856  V  1.75p.25p.75/(d  1.25g^)  (p^)  (RH  &  SH) 

(2) 

where: 

L  = 

ft ,  V  =  fi/sec ,  d  =  in ,  p  =  Ibm/ft^ ,  AP  =  psi  ,p  =  .0011886  ft2/sec , 

gc  = 

=  32.21bmfi/lbfsec2 

(d) 

m  =  pVA  mass  flow  rate 
m  =  .326pd2v  Ibm/min  (round  hole) 

(3a) 

m  =  .415d2v  Ibm/min  (square  hole) 

(3b) 

flow  rate  =  q  =  VA 

q  =  2.45d2V  (round  hole)  (gpm) 

(4a) 

q  =  3.12d2v  (square  hole)  (gpm) 

(4b) 

Convection  coefficient 

Nusset  number,  and  colbum  equation 

(a) 

Nu  =  .0223  Re-*  ft  -333  -  colbum  equatimi  where  assumptions 
(.7<ft<160,L/D>10) 

(b) 

Pr  =  v/a  (Prandtl  number) 

ft  =  Ratio  of  momentum  and  thermal  diffiisivities 

(c) 

a  =  K/pCp  (thomal  diffiisivi^  (Re  >  2300) 

(d) 

h  =  NuK/d  (convection  coefficient) 
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substituting  (a),  (b),  and  (c)  into  (d) 

h  =  .0223p-8V-8Cp-333K.667/(^.467d^)  (RH&SH)  (5) 

wha«: 

V  =  ft/sec ,  d  =  in .  p  =  IbnVft^ ,  p  =  .0011886  ft^/sec .  Cp  =  Btu/lbm  -T , 

K  =  Btu/hr.-ft  -T,  h  =  Btu/hr.- ft2  -T 

FILM  DROP  AND  FLUID  CHANGE  IN  TEMPERATURE 
Film  drop: 

(a)  Q  =  hAsATfd  (ATfd  =  Ts-  Tm)  (ten^peratuie  difference  between  the  fluid 
and  the  copper) 

(b)  Ag  =  ttdL  (surface  area  of  hole)  (round  hole) 
substituting  (b)  into  (a) 

ATfd  =  2.865(y(hd)  (round  hole)  (6a) 

(c)  Ag  =  4dL  (surface  area  of  hole)  (square  hole) 
substituting  (c)  into  (a) 

ATfd  =  2.25Q/(hd)  (square  hole)  (6b) 

where: 

h  =  Btu/hr.-  ft2  -T ,  d  =  in ,  Q  =  Btu/min ,  ATfd  =  *F 

Fluid  change  in  temperature: 

(a)  Q  =  rnCpATfi  (ATg  =  T  ouqxit  ■  "^input  fluid) 

ATfl  =  (y(.326pd2vq,)  (round  hole)  (7a) 

ATfl  =  (y(.415pd2VCp)  (square  hole)  (7b) 

where: 

V  =  ft/sec,  Cp  =  Btu/lbm  -*F ,  d  =  in ,  Q  =  Btu/min ,  ATfl=  T 


SOLUTION  FOR  THE  FINAL  SELECTED  GEOMETRY 

Geometry  of  the  cooling  hole  (round),  type  of  fluid  (Coolanol™  20),  and  the 
optimum  diameter  (d  =  .174")  have  been  established  based  on  thermal  performance  and 
practical  mechaitical  considerations.  Also,  it  has  been  established  that  the  fluid's  change 
in  temperature  (ATg  )  is  the  dmninating  force  in  the  thermal  analysis,  and  the  film  drop 
has  less  of  an  impact  The  above  equations  can  be  solved  in  terms  of  pressure  and  ATg, 
since  the  geometry,  diameter,  and  fluid  properties  (100  *F)  have  been  determined. 

Geometric  factors 

Acu  =  conductor  cross  section  =  (s-2t)2  -  icd2/4  =  .02376  (in^) 
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Acool  =  Jtti2/4  =  .0238  (in2) 

Pf  =  packing  factOT  -  Acu/At  =  .7604  -  12.56d2  =  .3801 

Heat  generation 

from  equation  (la) 

Q  =  1.941/(.047524  -  .785d2)  Btu/min  (8) 

substituting  d  =  .174"  into  equation  (1) 

Q  =  81.7  Btu/min 

Head  loss  and flow  rates 
from  equation  (2) 

AP  =  1.9856  V  1.75^.25p.75/(d  1.25gc)  (p^)  (9) 

substituting  d  =  .174",  p  =  55.2  Ibm/  ft^.  p  =  .0011886  Ibm/ft-sec  into  (2) 

AP  =  2.09V  1-75  (psi) 

from  equation  (3a) 

m  =  .326pd2v  Ibm/min  (10) 

substituting  d  =  .174",  p  =  55.2  Ibm/  into  (3) 
m  =  .545V  Ibm/min 
firom  equation  (4a) 

Q  =  2.45d2v  (gpm) 
substituting  d  =  .174"  into  (4a) 

Q  =  . 0742V  (gpm)  (11) 

Fluid  change  in  temperature: 
from  equation  (7a) 

ATfl  =  Q/(.326pd2vq,)  CF) 

substituting  d  =  .174",  p  =  55.2  Ibm/  ft^.  Cp  =  .47  Btu/lbm  -*F  into  (7a) 

ATfl=319.16/V(*F)  (12) 

The  first  of  two  assumptions  that  were  applied  in  order  to  simplify  the 
calculations,  was  that  the  fluid  properties  and  copper  resistivity  (taken  at  100  ^  would 
have  a  negligible  affect  on  ATg  as  the  temperature  changes  through  the  coil.  Hg.  C-1. 
summarizes  the  change  in  tenq)erature  of  the  fluid  (ATg)  varying  over  a  range  of  average 
temperatures.  Since  the  increased  copper  loses  at  higher  temperatures  (as  a  result  of 
higher  resistivity)  are  almost  equally  offset  by  lower  coolant  viscosity  (and  thus  higher 
flow  at  a  given  pressure  differential),  ATg  is  relativity  insensitive  to  fluid  temperature. 
The  assumption  is  valid.  The  second  assun^tion  involved  the  pressure  drop  across  the 
coil.  The  pressure  drop  was  calculated  to  be  a  straight  pipe  unaffected  by  the  coil's 
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curvature.  Fig.  C— 2.  displays  how  the  change  in  ten:q>erature  of  the  fluid  is  affected  by 
the  increased  pressure  drop  due  to  the  curvature  of  the  coil.  There  is  a  5  T  difference  in 
the  change  in  temperature  of  the  fluid  between  the  straight  pipe  assumption  and  the  actual 
curved  pipe  conduction. 

Eressure  drop  across  a  coil  (including  the  curvaturel 

AP  =  [.32/Re  *25  +  .048(d/D)-5]  (L/d)  (pV2/2gc)  (for  R©  =  1.5x10^)  (13) 

D  =  average  diameter  of  the  coil  =  8.51" 

substituting  D  =  8.51",  d  =  .174",  p  =  55.2  Ibm/  ft^.  p  =  .(X)l  1886  Ibm/ft-sec, 

L  =  80'  into  Eq.  5 

AP  =  [.3614V1*25  +  .226V2]  (14) 

where: 

L  =  ft ,  V=  ft/sec ,  d  =  in ,  p  =  Ibm/ft^ ,  AP  =  psi ,  p  =  .0011886  ft^/sec ,  D  =  in 
gc  =  32.2  Ibm  ft/lbfsec2 
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Fig.  C-1.  Fluid  AT  Vs.  Average  Ruid  Temperature 
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APPENDIX  D 

DETAILS  OF  STRUCTURAL  ANALYSIS. 

A  section  view  of  the  motor  is  shown  in  Fig.  D— 1.  The  inner  shaft  is  supported  by 
two  pairs  of  bearings  23.8  inches  apart  These  bearings  transmit  the  weight  of  the  mner 
shaft  to  the  outer  shaft  The  outer  shaft  is  also  supported  by  two  pairs  of  bearings  just 
outboard  of  the  inner  shaft  bearings.  The  outer  shaft  bearings  support  the  total  weight  of 
both  shafts  in  the  housing.  The  housing  is  assumed  to  be  held  at  the  bolted  flange 

connections  of  the  service  and  drive  end  plates. 

The  inner  shaft  material  is  1018  hot  roUed  steel.  The  outer  shaft  material  is  304 
stainless  steel,  with  iron  bars  and  copper  coils  adding  weight  but  not  structural  stiffness. 
The  housing  is  composed  of  316  stainless  steel,  A36  steel  and  A53  steel  pipe.  Table  D-1 
contains  the  Tnaterial  properties  for  the  components.  The  nomenclature  of  the  variables 
referred  to  is  alsf*  listed.  The  endurance  strength  for  each  material  has  been  calculated. 


TABLE  D-L  Material  Properties: 


Material 

Sv  (psi) 

E  (psi) 

5 

1018  Steel, 

Rod 

32,000 

58,000 

.3 

29E6 

29 

304 

Stainless 

Steel,  P*ipe 

30,000 

75,000 

.3 

28E6 

.29 

316 

Stainless 

Steel,  Pipe 

35,000 

85,000 

.3 

28E6 

.29 

A53  Steel, 

Pipe 

35,000 

60,000 

.3 

29E6 

.29 

A36  Steel, 

Plate 

36,000 

58,000 

.3 

29E6 

.29 

Iron,  Bar 

for  the  iron  and  copper  only 
the  density  was  used  in  the 
model 

.26 

.32 
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Fig.  D-1.  Section  View 


NOMENCLATURE: 

As  =  Shear  Area,  sq  in. 

At  =  Tensile  Area,  sq  in. 

D  =  Outer  Diameter,  in. 

p  =  Inner  Diameter,  in. 

F  =  Belt  Tension,  lb 

Fbn  =  Axial  Force  of  Bolt  due  to  Moment,  Mbn>  lb 
K  =  Endurance  Modification  Factors 

Kb  =  Bending  Stress  Concentration  Factor 

Kc  =  Combined  Stress  Concentratimi  Factor 

Kl  =  Load  Cotmection  Factor 

Kt  =  Torque  Stress  Concentration  Factcn* 

L  =  Distance  from  restraint  to  load,  in. 

M  =  Total  Moment,  in-lb 

Mbn  =  Moment  of  Bolt,  in-lb 

M'  =  Corrected  Moment,  in-lb 

N  =  Rotational  Speed,  rpm 

P  =  Power,  hp 

p  =  axial  force  on  bolt,  lb 

r  =  Radius,  in. 

Sa  =  Alternating  Strength,  psi 

Sc  =  Fully  Corrected  Endurance  Strength,  psi 


Se' 

=  Endurance  Strengdi,  psi 

Sm 

=  Mean  Strength,  psi 

Sse 

=  Shear  Endurance  Strength,  Psi 

Ssy 

=  Shear  Yield  Strength,  psi 

St 

=  Tensile  Strength,  psi 

Sy 

=  Yield  Strength,  psi 

s 

=  Shear  force  on  bolt,  lb 

T 

=  Torque,  in-lb 

wt 

=  weight,  lb 

zn 

=  Height  of  Bolt  to  center,  in. 

zn’ 

=  Corrected  Height,  in. 

a 

=  Angle,  degrees 

a 

=  Principle  Stress,  psi 

=  Poission's  Ratio 

c 

=  Von  Mises  Stress,  psi 

5 

=  Density,  Ib/cu  in. 

<ya 

=  Alternating  Stress,  psi 
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£  =  Summation 


Gin  =  Mean  Stress,  psi 
T  =  Shear  Stress,  psi 

ENDURANCE  STRENGTH  CALCULATIONS 
Endurance  Strength 


Sc 

=  Kl*K2*K3*Se* 

(REF.  11) 

K1 

=  Surface  Factor 

K2 

=  Size  Factor,  .869d-*097 

(select  largest  diameter) 

K3 

=  Reliability  Factor 

Inner  Shaft.  1018  Steel 

K1  =  Surface  Factor,  hot  rolled,  =  .7 

K2  =  Size  Factor,  .869  7.4-  097,  =  .72 


K3 

=  Reliability  Factor,  .999, 

=  .75 

Se' 

=  l/2St  =  29,000  psi  Se  = 

msiM 

Outer  Shaft.  304  Stainless  Steel 

K1 

=  Surface  Factor,  hot  rolled,  = 

.7 

K2 

=  Size  Factor,  .869  17-  097, 

=  .66 

K2 

=  Size  Factor,  .869  5.1-*097, 

=  .74  @  node  46 

K3 

=  Reliability  Factor,  .999,  = 

.75 

Se’ 

=  1/2  St  =  42,5000  psi 

Se  =  14.726  psi 

Sc  =  16.511  psi  @  node  46 

Housing.  316  Stainless  Steel 

K1 

=  Surface  Factor,  hot  rolled,  = 

.7 

K2 

=  Size  Factor,  .869  ll-*097. 

=  .69 

K3 

=  Reliability  Factor,  .999,  = 

.75 

Sc 

=  1/2  St  =  42,5000  psi 

Se  =  15.396  nsi 

Housing.  A53  Steel 

K1  =  Surface  Factor,  hot  rolled,  =  .7 
K2  =  Size  Factor,  .869  18-097,  =.66 
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K3  =  Reliability  Factor,  .999  =.75 


Sc’  =1/2  St  =  30,000  psi  Sc  =  10.395  psi 

Housing.  A36  Steel 

K1  =  Surface  Factor,  hot  rolled,  =  .7 

K2  =  Size  Factor,  .869  21.6'*09^»  =.65 

K3  =  Reliability  Factor,  .999  =.75 

Sc'  =1/2  St  =  29,000  psi  Sc  =9.896  psi 

The  two  shafts  and  the  housing  were  modeled  using  a  PC  based  finite  element 
package,  MSCNPAL2.  The  finite  element  model  comprised  of  circular  beams  represented 
the  geometry  of  the  Uiree  components  were  used.  The  flanged  connections  were  modeled 
as  rigid  coimections.  The  modeling  of  the  inner  shaft  with  circular  beams  is  very 
accurate.  However,  the  analogy  is  flawed  when  modeling  the  outer  shaft  and  the  housing 
for  pieces  such  as  the  stub  shaft  or  transition  flange.  Those  pieces  are  disks,  not 
cylinders.  To  properly  model  those  shapes  would  require  detailed  modeling  of  them  as  3- 
D  tvick  elements. 

The  shafts  were  subject  to  two  types  of  loads.  The  mean  loads  are  due  to  the 
torque  from  the  motor.  The  alternating  loads  are  fiom  the  shaft  bending  under  its  weight 
or  a  side  force.  The  sQess  results  were  separated  into  mean  and  alternating  cases.  Three 
load  cases  were  used; 

a)  The  steady  torque  from  the  300  hp  at  1200  rpm,  15,756  in-lb. 

b)  The  weight  of  the  components  plus  side  forces  (e.g.  from  belt 

force  on  outer  shaft) 

c)  Load  case  a  and  b  combined. 

The  bearings  were  modeled  as  very  stiff  ^rrings  connecting  the  components  in 
translations  but  not  in  rotation.  The  tCHque,  15,756  in-lb  due  to  300  hp  at  1200  rpm,  was 
placed  at  the  middle  of  the  iron  bars  between  the  bearings  on  each  shaft  The  rotations 
for  the  housing  and  the  shafts  were  fixed  along  their  centoline,  x-axis.  All  water  boxes, 
bearings,  carriers  and  pulleys  were  modeled  as  cylinders  with  estimated  weights  and 
dimensions.  See  Rgs.  D-2,  D-3,  and  D-4  for  individual  finite  element  (FE)  models 
showing  the  iimer  shaft  outer  shaft  arid  housing,  respectively. 
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LOCATION  5 
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Fig.  D-3.  Outer  Shaft 


Rg.  D-4.  Housing 
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Fcjr  the  dynamic  cases,  the  FE  model  was  held  from  translations  and  was  held 
from  rotating  in  all  directions  as  the  drive  end  so  as  to  have  a  stable  model.  A 
satisfactory  result  would  have  the  lowest  natural  frequency  well  above  the  maYimiim 
operating  speed. 

The  lowest  natural  frequency  of  the  assembly  was  at  69  Hz  or  4140  ipm,  which  is 
more  than  three  times  higher  than  the  maximum  operating  speed.  The  shafts  and  housing 
supported  as  described  should  be  acceptable. 

At  locations  where  changes  in  diameter  occur,  the  stresses  were  multiplied  by 
stress  concentrations  factms,  torsion  ot  bending,  respectively  for  the  first  two  cases.  For 
the  combined  case,  the  larger  of  the  two  factors  was  used.  To  calculate  the  actual 
combined  stress  concentration  factor,  Kc,  is  very  involved  and  assuming  the  larg^  of  the 
two  separate  cases  is  to  err  conservatively.  See  tables  D— 2  through  D-4  for  the  stress 
concentration  factors  applied  to  the  three  load  cases. 

TABLE  D-2.  Inner  Shaft  Stress  Concentration  Factors 


NODE 

D/d 

r/d 

Kt 

Kb 

Kc 

8 

1.07 

.01* 

1.8 

2.58 

2.58 

9 

1.10 

.01* 

1.8 

2.7 

2.7 

10 

1.35 

.06 

1.65 

1.9 

1.9 

11 

1.30 

.01* 

2.15 

2.8 

2.8 

12 

1.02 

0* 

1.8 

2.3 

2.3 

14 

1.11 

0 

1.9 

2.7 

2.7 

16 

1.29 

.01 

2.15 

2.8 

2.8 

17 

1.59 

^08 

1.55 

1.7 

1.7 

18 

1.10 

.01* 

1.8 

2.7 

2.7 

26 

1.07 

.01* 

1.8 

2.58 

2.58 

TABLE  D— 3.  Outer  Shaft  Stress  Concentration  Factors 


NODE 

D/d 

r/d 

Kt 

Kb 

Kc 

50 

1.04 

.01 

1.8 

2.6 

2.6 

36 

1.65 

.14 

1.3 

1.5 

1.5 

44 

1.77 

.19* 

1.25 

1.4 

1.4 

52 

1.05 

.01* 

1.8 

2.6 

2.6 
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TABLE  D-4.  Housing  Stress  Concentration  Factors 


NODE 

D/d 

r/d 

Kt 

Kb 

Kc 

62 

1.6 

.09 

1.45 

1.7 

1.7 

63 

1.2 

.06 

1.55 

1.9 

1.9 

66 

1.1 

.05 

1.3 

1.85 

1.85 

70 

1.1 

.05 

1.3 

1.85 

1.85 

71 

1.5 

.13 

1.4 

1.6 

1.6 

*  -  Radius  is  determined  by  bearing  company  specs, 
assumed  1/32"  radius 

These  values  were  calculated  from  diagrams  from  Ref.  1 1. 

The  results  from  the  frist  two  load  cases  generated  mean  Von  Mises  stresses  and 
alternating  Von  Mises  stresses,  respectively.  The  alternating  and  mean  stresses  were 
multiplied  by  the  appropriate  stress  concentration  factors  and  then  plotted  on  the  stress 
axes  of  a  Soderberg  diagram.  The  fatigue  safety  factor  was  calculated  by  dividing  the 
alternating  strength  by  the  altmiating  stress  (of  the  mean  strength  by  the  mean  stress). 
The  alternating  and  mean  strengths  are  draived  from  a  line  drawn  firom  the  origin  through 
the  plotted  stress  point  and  intersecting  the  Soderberg  line.  The  results  from  the  third 
load  case  generated  maximum  Von  Mises  stresses  fOT  the  static  load  condition.  These 
stresses,  after  being  multiplied  by  Kc,  were  compared  to  the  yield  sttength  for  the  static 
factor  of  safety  per  reference  1 1. 

It  was  assumed  that  direct  drive  coupling  would  be  used  to  transmit  power  from 
the  nsotors  inner  shaft  to  a  dynamometer  since  calculations  proved  that  the  iimer  shaft 
could  not  support  the  side  loads  imposed  by  a  belt  type  drive. 

To  transmit  power  from  tiie  motm*  te  a  dynamometer,  the  outer  shaft  was  assumed 
to  have  a  belt-pulley  type  connection.  For  tius  application  a  custom  pulley  will  be 
required.  A  Browning  type  pulley  (12"  diameter,  104  lb)  provided  a  weight  estimate.  A 
belt  force  of  3939  pounds  was  applied  at  the  drive  end  of  die  shaft 

The  stress  and  safety  factors  for  the  three  components  are  displayed  in  Table  D-5. 
The  deflections  are  in  Table  D-6.  The  Soderberg  graphs  for  the  three  components  are 
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shown  in  Figs.  D-5,  D-6,  and  D-7.  Hand  calculations  verifying  stress  levels  for  the 
inner  and  outer  shafts  are  also  presented. 

TABLE  D-5.  Stress  And  Safefy  Factors 


COMPONENT 

CASE 

NODE 

VONMISES 

STRESS 

K 

MOD. 

STRESS 

SAFETY 

FACTOR 

INNER  SHAFT 

a 

18 

5399 

1.8 

9718 

- 

b 

18 

580 

2.7 

1566 

2.4 

c 

18 

5430 

■ 

14661 

2.2 

OUTER  SHAFT 

a 

52 

1796 

1.8 

3233 

b 

52 

2810 

2.6 

7306 

1.9 

c 

52 

3131 

2.6 

8141 

3.7 

HOUSING 

a 

• 

b 

70 

182 

1.9 

337 

46 

c 

70 

182 

1.9 

337 

104 

TABLE  D-6.  Deflections 


COMPONENT 

NODE 

MAX.DEFL 

NODE 

DEFLECT.  BETWEEN 

BEARINGS 

INNER SHAFT 

1 

l.lE-3 

17 

1.6E-4 

OUTER  SHAFT 

49 

3.6E-3 

44 

6.4E-5 

HOUSING 

72 

2.0E-4 

68 

8.2E-6 
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PRELIMINARY  CALCULATION  OF  HOMOPOLAR  MOTOR  INNER  SHAFT 
Assumptions: 

Power  to  be  transmitted  is  300  hp. 

Shaft  material  is  1018,  hot  rolled  steel 

A  pulley  was  found  from  the  Browning  Catalog  to  transmit  300  hp  on  a  2  3/4"  diameter 
shaft,  a  Poly  V  drive  with  15"  diameter.  The  weight  of  the  pulley  and  half  the  belt  is  1 13 
lb. 

Calculations: 

T  =  63026  P  /  N  =  63025*300/1200  =  15,756  in-lb 
F  =  2  T  Kl  =  L5  for  pulley  connection 
=  2*15756*1.5/15  =  113  =  3264  lb 
Reverse  Bending 

this  calculation  is  based  on  belt  tension  acting  at  the  right  end  of  shaft,  the  longest 
unsuppcH^  span.  Hie  calculation  is  simplified  by  assuming  a  constant  diameter  of  2 
3/4"  from  right  bearing  to  end  of  shaft 


M  =  F  L  =  3264*12  =  39,168  in-lb 

sa  =  32  M  d4=  19.190 psi  (REFll) 
nCD^-d^) 

Se=  10,962  psi 

Sa  >  Se  therefrne  shaft  is  too  small 
for  ovohung  load 

STRESS  CALCULATION  VERIHCATION 
Assumptions: 

Power  to  be  transmitted  is  300  hp. 

The  weight  of  half  the  coupling  is  15  lb. 
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The  weight  of  the  pulley  and  half  the  belt  is  102  lb. 

Shaft  is  cantilever  beam,  double  bearings  act  as  fixed  end  condition,  with  uniform  load 
due  to  weight  and  a  concentrated  load  due  to  weight  of  pulley  or  coupling. 

INNER  SHAFT 

Calculations: 

I  =  JC  (  d4  -  d^)  /64  =  n  (2.9534-.3754)  /64  =  3.7  in^ 

J  =  jc  (  d4  -  d4)  /32  =  jc  (2.9534-.3754)  /32  =  7.5  in^ 
w  =  d  Jt  (  d2  -  d2)  /4  =  .28  JC  (2.953^-375^)  lA  =  1.9  Ib^ 

Steady  Stress 

t  m  =  T  r  /  J  =  15756(1.48)  /7.5  =  3120psi 

3117  psi  fix)mMSCyPAL2 


Altemating  Stress 


I 


151b 

i 


M  =  FI 
=  15(13.4) 
=  201  in-lb 


M  =  w  ^  /  2 
=  1.9  (13.42) /2 
=  171  in-Ib 


s  =  M  c  /  I  =  (201+171)1.48/  3.7  =  147  psi 


Maximum  Von  Mises  Stress 

s'  =  (s2  +3tm2)2  =  (1472  +  3(31202))2 

s' =  5406  Psi 

5400  psi  finm  MSC/PAL2 
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OUTER  SHAFT 

Calculations: 

I  =  It  (  d4  -  d4)  /64  =  1C  (5.124-4.134)  /64  =  19.5  in4 
J  =  ic  (D4-d4)/32  =  ic  (5.124^.134) /32  =  38.9  in4 
w  =  d  It  (  d2  -  d2)  /4  =  .28  It  (5.122-4.132)  /4  =  2.0  ll/m 
Steady  Stress 

t  m  =  T  r  /  J  =  15756(2.56)  /38.9  =1037  psi 

1037  psi  fiom  MSCyPAL2 


Alternating  Stress 


I 


151b 

i 


M  =  FI 
=  102(4.9) 

=  499.8  in-lb 


M  =  / 2 

=  2(7.4^)/2 
=  54.8  in-lb 


s  =  M  c  /  I  =  (500+55)2.56/ 19.5  =  72.8  psi 


Maximum  Von  Mises  Stress 

s*  =  (s2  +3tm^)2  =  (732  +  3(10372))  2 

s‘=1798  psi 

1796  psi  fiom  MSCyPAL2 
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BOLT  LOADING,  STRESSES  AND  SAFETY  FACTORS 

The  bolt  loads  for  the  major  connections  of  the  outer  shaft  and  the  housing  were 
calculated  for  shear  and  axial  stresses.  All  the  bolts  arc  assumed  to  be  pieloaded  so  that 
none  will  be  subject  to  cyclic  loading.  A  high  preload  improves  the  fatigue  resistance  of 
a  bolted  coimection.  The  bolt  coimection  locations  are  shown  in  Hgs  D-3  and  D-4. 

Note,  this  model  was  developed  to  determine  the  stress  and  deflections  of  large 
components,  which  was  accomplished  by  using  coarse  circular  beam  elements.  It  was  not 
intended  to  be  used  to  solve  areas  of  local  stress  on  a  small  scale.  These  results  for  the 
bolts  are  only  useful  as  a  rmigh  guide  for  sizing  and  numbers  required.  See  tables  D-7 
through  D-12  for  the  results  of  the  calculations  of  the  bolt  loading,  stresses  and  safety 
factors. 


Shear  Stress 
t_(»+T/d) 

As 

Ssy  =  .577*Sy  (REF.  11) 

tmax  =  ((s/2)2  + 

smax  =  tmax  •¥  s/2 

Assume  all  bolts  are  preloaded  so  that  there  will  be  no  cyclic  loading 
on  any  bolts. 


Axial  Stress 
At 


Neutral 

Axis 
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Calculation  of  Axial  Bolt  Load from  a  Moment 

The  farce  on  each  bolt  is  proportional  to  the  distance  ftom  the  neutral  axis.  It  is 
assumed  that  the  ration  of  force  on  each  bolt  was  a  quadratic  of  the  Hktancr 

For  an  8  bolt  connection: 

zi=r  22  =  rsina  Xz  =  2*2i' +  4*22 

^  ~  I  *1^  =1  Z2'  =  .5,  for  any  radius 


ThcrefOTcXz  =  2*l  +  4*.5  =  4 

M’  =  M/l2  Mbn  =  M’  ♦  2n'  Fbn  =  Mbn  /  Zn 

In  table  D-7  axial  bolt  load  was  only  calculated  for  the  bolt  farthest  away  ftom 
the  center,  Fbi,  the  worst  case;  Fbi  =  Mbl  /  z\ 

For  any  insulated  bolts,  it  is  assumed  that  only  die  shank  is  Mmed  down  and 
insulated  and  that  the  shank  diameter  is  not  less  than  the  minor  diameter. 

BOLT  LOADING,  STRESSES  AND  SAFETY  FACTORS  OUTER  SHAFT 

Group  1 

Assume  1/4-20  UNC  bolts,  8  per  connection 

Sy  =  55,000  psi  Ssy  =  31,735  psi  Grade  2  Steel 

As  =  .027sqin  At  =  .032sqin 
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TABLE  D-7.  Total  Forces  At  Connection 


Bolt 

Axial 

Shear 

Torsion 

Moment 

P 

s 

T 

M 

Location 

LB 

LB 

IN-LB 

IN-LB 

1 

0 

65 

0 

177 

2 

0 

904 

0 

2170 

3 

0 

962 

0 

6181 

4 

0 

1140 

0 

8491 

5 

0 

0 

15760 

0 

TABLE  D-8.  Stresses  Per  Bolt 


AXIAL 

Fb 

t 

PSI 

s 

PSI 

tmax 

PSI 

smax 

PSI 

FS 

t 

FS 

s 

1 

22 

301 

691 

458 

804 

69 

68 

2 

103 

4185 

3229 

4486 

6100 

7 

9 

3 

252 

4454 

3229 

5948 

9890 

5 

6 

4 

260 

5278 

8139 

6665 

10734 

5 

5 

5 

0 

9728 

0 

9728 

9728 

3 

6 

Group  2 


Sy  =  170,000  psi  Ssy  =  98,090  psi  Unbrako  Series  1960 
TABLE  D-9.  Total  Forces  At  Connection 


BOLT 

AXIAL 

SHEAR 

TORSION 

MOMENT 

P 

s 

T 

M 

LOCATION 

LB 

LB 

IN-LB 

IN-LB 

6 

0 

1356 

15760 

21100 

7 

0 

1451 

15760 

23580 

8 

0 

1494 

15760 

23590 
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TABLE  D-10.  Stresses  Per  Bolt 


BOLT 

AXIAL 

Fb 

t 

PSI 

s 

PSI 

tmax 

PSI 

smax 

PSI 

FS 

t 

FS 

s 

6 

15258 

20288 

18322 

28467 

5 

6 

7 

813 

16781 

25409 

21048 

33753 

5 

6 

8 

1311 

23131 

40955 

30893 

51370 

3 

3 

BOLT  LOADING,  STRESSES  AND  SAFETY  FACTORS  HOUSING 

Assume  1/4-20  UNC  bolts,  8  per  coimecticm 

Sy  =  55,000  psi  Ssy  =  31,735  psi  Grad  2  Steel 

As  =  .027  sq  in  At  -  .032  sq  in 

TABLE  D-11.  Total  Forces  At  Connection 


BOLT 

AXIAL 

SHEAR 

TORSION 

MOMENT 

P 

s 

T 

M 

LOCATION 

LB 

LB 

IN-LB 

IN-LB 

1 

0 

392 

0 

399 

2 

0 

1132 

0 

4938 

3 

0 

1232 

0 

21290 

TABLE  D-12.  Stresses  FCr  Bolt 


AXIAL 

Fb 

t 

PSI 

s 

PSI 

tmax 

PSI 

smax 

PSI 

FS 

t 

FS 

s 

1 

10 

1815 

300 

18322 

1971 

17 

28 

2 

119 

5241 

3718 

21048 

7420 

6 

7 

3 

374 

5704 

11672 

30893 

13996 

4 

4 
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APPENDIX  E 
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N1  IKILt,  /U  UUKUHfe.lkK. 

COKPOUf^O  NO  N674-70,  '^ARKER 


nUKMAL  l^UnLTUL  I  JVC,  HUMUrULAK, 

CONTRAfiOTAJlNG,  1^5  HP  MOTOR 


ASSEMBLY 

SECTION 


T 


di 


< 


127-19392-6 


INNER  ROTOR 
ASSEMBLY 
SECTION 


o 


< 


u 


< 


SERVICE  END 
OUTER  ROTOR 
RETURN  RING 
DETAILS 


0 


< 


DRIVE  END 
STUB  SHAFT 
DETAILS 


i27H9992-^ 


BENDS  DETAILS 


< 


CONDUCTOR 

DETAILS 


BEAR  ING 

retainer 

DETAILS 
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